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Abstract 

A density-matrix formalism is developed based on the one-particle density-matrix of a single- 
determinantal reference-state. Unlike traditional density-functional-theory approaches, the v- 
representable problem does not appear in the proposed method, nor the need to introduce function- 
als defined by a constrained search; furthermore, the necessary conditions for a one-particle density 
matrix to come from a single determinant are known, so they can be expressed as constraints when 
minimizing variational-energy functionals. The correlation-energy functionals are not universal, in 
the sense that they depend on the external potential. Nevertheless, model systems can still be 
used to derive universal energy-functionals. For example, the Colle-Salvetti functional is shown 
to be compatible with the proposed method. In addition, the correlation-energy functionals can 
be partitioned into individual terms that are - to a varying degree - universal; yielding, for exam- 
ple, an electron gas approximation, where the gas in not assumed to be uniform. Variational and 
non-variational energy functionals are introduced that yield the target state when the reference 
state - or its corresponding one-particle density matrix - is constructed from Brueckner orbitals. 
Using many-body perturbation theory, diagrammatic expansions are given for the non-variational 
energy-functionals, where the individual diagrams explicitly depend on the one-particle density- 
matrix. Non-variational energy-functionals yield generalized Hartree-Fock equations involving a 
non-local correlation-potential and the Hartree-Fock exchange; these equations are obtained by 
imposing the Brillouin-Brueckner condition. The same equations - for the most part - are ob- 
tained from variational energy-functionals using functional minimizations, yielding the (kernel of) 
correlation potential as the functional derivative of correlation-energy functionals. Approximations 
for the correlation-energy functions are introduced, including a one-particle-density-matrix vari- 
ant of the local-density approximation (LDA), a variant of the Lee- Yang-Parr (LYP) functional, 
and a variant of the correlation-energy functional used in the so-called B3LYP approach. A brief 
comparison with the exact SCF theory by Lowdin is presented. 
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I. INTRODUCTION 



Many variants of density functional theory (DFT) [l|, y, |3|, y , |5|, |6[ share common features 
with the Hartree-Fock approach. In particular, the Kohn-Sham method Q uses orbital 
equations that appear similar, in certain respects, with the ones from Hartree-Fock theory. 
On the other hand, unlike the Hartree-Fock wavefunction, the Kohn-Sham determinantal 
state shares only a common density with the exact wavefunction, and is not considered an 
approximation of the ground-state. Furthermore, in contrast to the exchange potential from 
Hartree-Fock theory, the Kohn-Sham exchange-correlation potential is local. 

DFT approaches that use hybrid functionals 0,0,0,^1 introduce a component of exact 
exchange-energy , w here justification, in part, for this modification comes from the adiabatic 
connection Iisl ll^ . [l5 . 16^ . yielding an approach that, again, has more similarities 
with Hartree-Fock, especially since its energy functional yields a non-local potential - the 
Hartree-Fock exchange - that depends on the one-particle density matrix of the Kohn-Sham 
determinantal- state. The most celebrated hybrid functional, B3LYP |^ , contains three 
parameters, two correlation-energy functionals, the Dirac-exchange functional (with a cor- 
rection), and, of course, exact exchange. The LYP density functional - a key component 
of B3LYP - is derived from the Colle-Salvetti correlation-energy functional where this 
functionals depends on a one-particle density-matrix, say 7, where 7 is from the Hartree- 
Fock reference state, and not an exact eigenstate, indicating a further evolutionary step of 
DFT methods towards a Hartree-Fock generalization with inclusion of electron correlation. 

In the Hartree-Fock Kohn-Sham approach [2^, the exchange energy is treated in an 
exact manner and the non-local, Hartree-Fock exchange-potential appear in the orbital 
equations. A generalization of this approach by Lindgren and Salomonson ^ yields, in 
addition, a nonlocal correlation potential and orbitals that, they believe, are very similar to 
Brueckner orbitals. Other workers also suggest that Brueckner and Kohn-Sham orbitals are 
very similar 



Brueckner orbital theory 
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24, H H 22, H Q, H 0, S Q, Q Q Q is 

a generalization of Hartree-Fock theory that utilizes a sing le-determinantal state that has 



38, 



39l |. Below we use this formalism 



the maximum overlap with an exact eigenfunction 
to develop a density-matrix theory, in which a variety of variational and non-variational 
energy-functionals are introduced that depend on the one-particle density-matrix, say 7. 



Jnlike othe r ap proaches, where 7 is the one-particle density-matrix of an exact eigenfunction 
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the introduced method - called reference- state one-particle density- 



matrix theory - has 7 arising from a single-determinantal reference state, where the energy- 
functionals yield the exact energy when 7 is the one from the Brueckner reference-state. 

One advantage that this one-particle, density-matrix approach has over traditional 
density-functional formalisms - or one-particle density-matrix formalisms - is that there 



is no w-representable pro 
a constrained search l43l. 



3leinjl|, 3, Q 



nor the need to introduce functionals defined by 
471 148| . Furthermore, the necessary conditions for a one-particle 
density matrix 7 to come from a single determinant are known, and they can be expressed 
as constraints when minimizing energy functionals that depend on 7 j2 1^ • 

Below, generalized Hartree-Fock equations are obtained containing the exact exchange- 
potential and a nonlocal correlation-potential, where these equations are obtained using 
the Brillouin-Brueckner condition - using non-variational energy functionals - and func- 
tional minimization - using variational functionals. Both variational and non-variational 
approached lead to the same correlation potential and generalized Fock-operator. (Cor- 
relation potentials from either approach are the same, in the sense that the (occupied) 
Brueckner-orbitals obtained from the variational correlation-potential differs only from the 



orbitals obtained from the non-variational one by unitary transformation. 
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56|, diagram- 



Using time-independent many-body perturbation theory 
matic expansions are given for the non-variational energy-functionals that are expressed in 
terms of orbitals and orbital energies. When restrictions are placed on the orbital energies, 
the individual diagrams are shown to explicitly depend on the one-particle density-matrix 
of the reference state. The diagrammatic expansions for the variational-energy functionals 
are presented elsewhere js^]. 

Kohn-Sham variants of DFT employ a universal exchange- correlation functional, inde- 
pendent of the external potential; approximations can be derived from model systems, where, 
in the vicinity of the model systems, the general form of the exchange-correlation functional 
is known. In contrasts, the correlation-energy functionals introduced below depend on the 
external potential, and are, therefore, in this sense, not universal. Nevertheless, as shown 
below, the correlation-energy functionals can be partitioned into individual terms that are 
- to a varying degree - universal; approximations can be derived from model systems. For 
example, the electron-gas correlation-energy can be used in an electron gas approximation, 



where, unlike the local density approximation (LDA) [7|, the gas in not assumed to be 
uniform. 

In addition, even without partitioning of the correlation-energy functionals - with the 
external-potential dependence intact - approximate functionals can still be derived from 
model systems. For example, as discussed below, the Colle-Salvetti functional [3| ^ derived 
from the helium atom - is a valid approximation within the proposed method. 



II. OVERVIEW 

Sec. IVl introduces four trial wavefunctions - say I^E'^'), where rj = I, II, III, and IV - that 
are defined with respect to a single-determinantal reference-state, say |$). The first trial- 
wavefunction is simply the target state of interest, say |\E'), with the single excitations 
removed. The second trial- wavefunction is defined with respect to the target state 

expressed by an exponential ansatz: (|\E') = e^\^)), where is generated by removing 

the single-excitation amplitudes .Si from the cluster-operator S. All of the trial states 

- including the third and fourth ones defined below - share the property that they contain 
no single excitations; furthermore, they generate the target state of interest |^) when their 
reference states satisfies (1$) = |0)), where |B) is the determinantal state constructed from 
occupied Bruckner orbitals. In other words, we have (I^E^q') = |^))- 

Using intermediate normalization, the exact energy of interest, say S, is just 
and it can be partitioned into two terms: the first-order energy, say - given by 

- and the correlation energy, say £^co['^'] ^ which, unlike S, is a functional of the reference 
state 1$); these concepts are briefly review in Sec. IIII Al 

Four non- variational energy-functionals, say i^rji"^]; are defined in an analogous way as S: 
{E^[^] = where they yield the exact energy for the Brueckner reference-state: 

{S = E^ilQ]), and these functionals can also be partitioned into two terms: = -Eii^] + 

-^co^[^])' where is the correlation- energy functional, and we have (£^co[0] = -E'c?[0]); 

the details are presented in Sec. El 

Because of the one-to-one correspondence between the set of determinant states, say 
{!$)}, and the one-particle density-matrices 0, 5^, say {7}, the correlation energy, say 
^^cob], and the correlation-energy functionals, say -Eco [7] ^ other functions and func- 

tionals of 7 - can be transfered into ones that depend on the one-particle density-matrix 7, 



as discussed in Sec. lVIIl"^ the correlation-energy and exact energy satisfy: (i^coH = E^o[t]) 
and {S = Ejj[t]), where r is the one-particle density matrix of the Brueckner reference-state 
16). 

Unlike the first-two trial wavef unctions, l^&ij') and 1^1^*^"'), that are generated by re- 
moving single-excitation amplitudes from the target state the third and fourth trial 
wayefunctio ns^ J^i:"') .^^4 1 ^7^. ^^' obtained by solving the _ coupled-cluster _ the or-^ 
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equations, respectively, in an approximate way - by neglecting the single-excitation portions. 

The coupled cluster formalism is briefly reviewed in Sec. IIII Bl A transparent perturbative 
treatment of the coupled cluster theory is presented in Sec. 1111 CI that is useful to obtain 
a perturbative expansion for the third correlation-energy functional E^™\ (This approach 
is compared to Lindgren's variant of the link diagram theorem in Sec. IIII EH Sec. IIII Dl 
reviews Rayleigh-Schrodinger perturbation theory that can be used to generate the fourth 
correlation-energy functional E^^J^ when the linked diagram theorem is not invoked and the 
single-excitation subspace is neglected. 

In Sees. IVIII Bl and IVIII diagrammatic expansions using many-body perturbation the- 
ory are presented for the correlation energy £^co[7] and the correlation-energy functional 
E^l^l'y]. The individual diagrams depend on the orbitals - both occupied and unoccupied 
- and the orbital energies - defined by the zeroth-order Hamiltonian. By using degenerate 
sets of occupied and unoccupied orbitals, and additional methods, it is demonstrated that 
diagrams can be defined that explicitly depend 7. 

In Sec. llXl an approach based on many-body perturbation theory is introduced, where the 
perturbation is partitioned into terms that depend on the external potential, say v, and the 
remaining portion that is v independent; the perturbation expansions for the correlation en- 
ergy Sco and the correlation-energy functionals E^l^ mirrors this partitioning, yielding terms 
that depend on v and, the remainder, called the electron gas terms, that are v independent. 
The electron gas terms are the only terms that contribute to Sco and E^l^ for an electron gas 
and, in most other the dominant portions of Sco and E^l^; they are also universal 

functions - independent of the external potential. Atomic, diatomic and molecular terms 
are defined in an analogous way are are obtained by further partitioning the perturbation 
into potential terms from the individual nuclei and selectively partitioning the perturbation 
expansions for Sco and E^^K 



Using the electron gas terms, an electron gas approximation is proposed in Sec |3 that 
is an ahernative to the LDA. Additional approximations are also considered in this Sec. 
including one that leads to the Colle-Salvetti functional [igj ]. 

Generalized Hartree-Fock equations are defined in Sec. IVlll where the exact Fock- 
operators, say generate the Brueckner orbitals, and these operators are defined by 

the trial wavefunctions \^%^) and the Brillouin-Brueckner condition, which is reviewed in 
Sec. IIVI In addition - from the one-to-one correspondence mentioned above - we can also 
write J^^''^; where, it is demonstrated that these operators are independent of r^; so, in ad- 
dition, we can omit the rj subscript and write J-'r- Solving the generalized Hartree-Fock 
equations permit the determination of the Brueckner orbitals, and the one-particle density- 
matrix, r, that is defined by these orbitals, permitting the determination of £ and i^coM, 
since, as mentioned above, they are given by Er^[T\ and i?^2'[r], respectively; A correlation 
potential, say f^^, is also defined, and satisfies the following identity: (JF^ = Ft + vl^), where 
Fr is the Fock operator, determined by r; it is constructed from the Brueckner orbitals. 

A variational formalism is presented in Sec. IXIf where energy functionals, say -^^[7], 
are defined using the same trial wavefunctions as in the non- variational case: (-E'r,[7] = 
{^^^^\H\^^^^) [(\E'*^'''|^^j;'')] ^), where the exact energy is generated for the Brueckner-state 
one-particle density- matrix: {£ = £'^[r]). These functionals are minimized subject to the 
constraint that the one-particle density-matrix comes from a single-determinantal state 
The functional derivative of Er)['~f] - with respect to the one-particle density-matrix 
- generates Ct''K^15^2), where these two-body functions are the kernels of generalized Fock 
operators, say and it is demonstrated that these operators are independent of 77; in ad- 
dition, these operators are - in a the sense mentioned in the introduction. Sec. |l]- equivalent 
to the non-variational operators, A correlation potential, say i)^^, is also obtained that 
is defined by its kernels which is given by the functional derivative of variational correlation- 
energy functionals, say E^l^[y]. The electron-gas and Colle-Salvetti functionals, mentioned 
above, are valid within the variational approach, permitting the determination of approxi- 
mate correlation potentials 0^^ by functional differentiation. 

A brief comparison with the Brueckner-orbital, exact SCF theory by Lowdin 2^ and 



Kobe 



38l | is presented in Appendix 1X1 



III. PERTURBATION AND COUPLED CLUSTER THEORY 



A. The exact and correlation energies 

We seek solutions of the time-independent Schrodinger equation, 

H\^)=S\^), (1) 
where is an eigenstate of the Hamihonian operator, 

H = '^[i\h\j]alaj + - "^[ijlk^alalaiaj, (2) 

ij ijkl 

and the integrals are written using chemist's notation 



mj] = m-h'^'M + mj], (3) 

[ij\kl] = J2 J V'r(xi)V^i(xi)rf2Vfc(x2)V^Kx2) dridr2, (4) 

where the spatial and spin coordinates, r and uj, are denoted collectively by x. 
The wavefunction of interest can be generated by a wave operator 

Q^\<l>) = 1^), (5) 

where |$) is a determinantal reference-state. 

The reference state |$) is completely defined by its occupied orbitals; we denote these 
orbitals by {i^o^ ^}', the set of unoccupied orbitals - the virtual orbital set - is denoted by 
{ipu^ The virtual set also determines the occupied set, since the two sets are 

orthogonal, and the union of the two sets is a complete set. Hence, {■0^^ ^} also determines 
1$). Unless stated otherwise, two sets of either occupied or unoccupied orbitals that differ 
by a unitary transformation are considered equivalent. 

We use the following orbital convention: Arbitrary orbitals are denoted by i and j; 
occupied orbitals are denoted by w, x, and y; virtual orbitals are denoted by r, s, and t: 

tpw.ipx.iJy e {V'o ^ $}, (6a) 

i^r^i^s^^t e {^pu $}, (6b) 

i)j,i)k e {ipo $} u {ipu (6c) 



Explicitly, our spin-orb itals ipi{'x) have the following form: 

Ai^) = Xiair)a{u), a = aoT /3, (7) 

where the spin and spatial portions are given by Xjo-(r) and cr(u;), respectively, and the spatial 
functions Xjo-(i") are permitted to be unrestricted - two spin orbitals do not, in general, share 
the same spatial function. 

In principle, |$) can be any determinantal state that overlaps with the reference state: 
7^ 0). However, our interest is often in cases where the target state is a ground 
state that is well described by a closed-shell reference-state |$), and the Hamiltonian is 
spin-free - it contains no spin coordinates. In these cases, instead of Eq. we often use a 
spatially restricted set of orbitals, given by 

ipiM = Xiir)(T{uj), a = a,(3, (8) 

so that 1$) is determined by a set of doubly-occupied spatial orbitals, denoted by {xo^ 
where this set also determines the virtual set {Xu^ Two sets of orbitals that differ by a 
unitary transformation are considered, again, as equivalent. 

By multiplying the Schrodinger Eq. (P) from the left by ($|, and requiring intermediate 
normalization to be satisfied, 

($1^) = = 1, (9) 

we get 

8={^\H\^) = E,[^]+S,om, (10) 

where the first-order energy is 

^i[$] = = Yl [w\{-'^V^) + v + l{J^-K^)\w], (11) 

and the Coulomb J$(r) and exchange -R'$(x) operators have their usual forms: 

[^I-^i-Ij] = Y [xx\ij], (12) 
[z\K^\j] = Yl (13) 
furthermore, the correlation energy £^co['^'], given by 

£eo[$] = (^I^^I^Q,), (14) 



is obtained from the correlation function: 



(15) 



where the orthogonal-space projector is 



g$ = i-|$)($|. 



(16) 



The first-order energy can also be written as 



(17) 



where the cl subscript indicates the closed portion - the fully contracted terms that, dia 



56, 



Appendix IHl presents 



grammatically speaking, have no external free-lines j5f 
partitioning or second-quantized operators into closed and open portions in a slightly differ- 
ent manner than is done by other authors. 

Similar to the first-order energy, for the correlation energy we have 



(18) 



where the correlation operator defined by 



X$ — — 1, 



(19) 



generates the correlation function | "^q^ when operating on the reference state: 



X<i.|$) = I^qJ. 



(20) 



The sum of Eqs. ffTTj) and (fTH|) gives the exact energy: 



cl • 



(21) 



We also write down the expression for the exchange-correlation energy: 



(22) 



where the exchange energy [$] is the last term on the right side of Eq. (Illj) : 



W\ 



(23) 



B. The linked cluster theorem 



The wave operator 17$ can be expressed in an exponential form [56|, |58|, |59|, |60|, |6l|, |62 , 



63|,|64, 



65| 
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1 



n^ = e'^ = i + s^ + -Si + -Si + 



(24a) 



where the cluster operator 5"$ can be written as a sum of one-, two-, and higher-body terms. 



and these amplitudes are defined by the following relations: 
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rwsxty 



(24b) 

(25a) 
(25b) 
(25c) 



which use the orbital convention given by Eqs. ® ; furthermore, the coefficients are required 
to satisfy exchange symmetry: 



3* 

'rwsxty 



Q C 

sxrwty '^rwtysxi 



(26a) 
(26b) 



The cluster operator 5"$ and its amplitudes S^ are invariant to unitary transformations of 
the occupied or virtual orbitals 

Since Si^ - given by Eqs. ()24b|) and - is open, only connected (cn) portions contribute 
to the correlation- and exact-energies, £^co['^'] and given by Eqs. (fTHj) and (PT|) . Therefore, 
we can write 



(27) 
(28) 



where the additional cn subscripts indicate that only the connected portions contribute 
contractions in which all 5"$ amplitudes are connected together by H. 



Eq. (^Hjl indicates that the closed part of (i70$)cn gives the energy of interest, £\ t 



ne open 



part is the mathematical statement of the linked-cluster theorem j56 
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63 



{Hn^)_^„ = 0. 



.op.cn - (29) 

Using the time-independent form of Wick's theorem the operator product 

HQij, can be written as a sum of zero-, one-, two- and higher-body excitations: 



Hn^ = £ + (Hn^)^ + (Hn^)^ + (nn^)^ + 



(30) 



where the notation (. . .)„ indicates that the n-body term within the brackets (. . .) is normal- 
ordered with respect to the |$) vacuum state. Substituting this expression into Eq. 
we get 



^)n\op,cn ~ 0- 



(31) 



n=l 



Since each term is linearly independent, the solution is 



KHn^)] =0, n^O. 



(32) 



This relation can be used to obtain the coupled-cluster equations, satisfied by the 5"$ am- 
plitudes, 5** 



Perturbation treatment of the linked cluster theorem 



For Rayleigh-Schrodinger perturbation theory 
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the 



wave- and cluster-operators are given by order-by-order expansions: 



Si 



0(1) , 0(2) , q{3) . 

"T -|- -|- 



4^) + ..., 



(33a) 
(33b) 



where 



n 



(0) 



(33c) 



Substituting these relations into Eq. ()24a|) and equating each order, we get the following 
identities jsy, : 

ni'^ = (34a) 

.(2) _ ^(2) , 1 /cW^^ 



^);' = + ^ [Si'>) , (34b) 
= sS' + Si^sS' + l-^[si^y, (34c) 



where the order from each term is defined as a sum of the superscripts, e.g., S'^^'S'^^' is a 
third order term. 

For a perturbative treatment, we partition the Hamiltonian into a zeroth-order Hamilto- 
nian Hq and a perturbation V: 

H = Ho + V, (35) 

where we require the reference state |$) to be an eigenfunction of Hq, and a one-body 
operator: 

i/o|<f) = Eo\<^), (36) 
Ho = J^^Q^aJaj. (37) 

The above zeroth-order Hamiltonian is defied by its matrix elements. We choose them 
by requiring the following relation to be satisfied: 

eij = eji = e* , (38a) 

where 

C = 0, (38b) 

c = (38c) 
4 = (Aia^), (38d) 

and the one-body operators, /* and /*, are determined by the reference state |$), but 
the dependence of /* and /* upon |$) is at our disposal; the orbital subspaces are, again, 
defined by Eqs. (jH)). 



Using the above choice, our zeroth-order Hamiltonian becomes 

^0* = Yl Cx^lax + Yl (39) 

where the appended $ superscript indicates that Hq now depends on the reference state |$). 

A hnked diagram expansion for Qs, and £^c.of3*Lis known to exist for a zeroth-order Hamil- 
tonian that is a diagonal, one-body, operator yHHQHHQQQ- 

A diagonal 

form for our one-body operator, Hq, is obtained when we choose its orbital sets - {ipo^ ^} 
and {i^u^ ^} ^ to satisfy the following conditions: 

mf!m = S^^C, (40a) 
(Aia^s) = Srset, (40b) 

where we denote these particular sets of orbitals by {i^o^ ^, /*} and {ipu^ /* }, indicating 
that they are uniquely determined by |$) and their one-particle operator, /* or /*. 
Using these orbitals, if* can be written as 

Hq = Y C(^l(^w + Y er4«r-, (41) 

where our partitioning is 

H = H^ + V^. (42) 

Since Hq is a second-quantized operator with no projection operators, it is appropriate 
to substitute the above partitioning into Eq. (j^ . yielding 

- W^*)op,cn=(^*^*)op,cn, (43) 

Substituting Eqs. (jH^ into PHj) . and equating each order, we have 



H^Sr =iV^n'r' , (44) 

/ op,cn V / op,cn 

where V$ contributes unity to the overall order of the rhs. (In a sense, we have: V^ = V^^\) 
Solving Eqs. ()44|1 and (jH^ in an iterative and sequential manner generates a linked 
diagram expansion for and a connected expansion for 5$, which is identical to the 

expansions obtained by Lindgen y, Q- 

The procedure begins by using Eq. (jHScj), 
which gives S^^ and fi^^'* from Eqs. pij) and (jMaj). Substituting ^l^^^ back into Eq. pij) 

^2) (2) 

gives Si , permitting the calculation of ^l)f, from ()34b|) . and so on. However, in order to 



get the same diagrams as Lindgren, the factorization theorem [52, |53|, |56|, [7^ l76|, l77|, [73, l7£ 
must by exploited to unfactorize the disconnected products that appear on the right side on 
Eq. (|34p . The diagrammatic expansions for the correlations and exact energies, £^co['^'] and 

are obtained from Eqs. (j27j) and (j28p . It is worth mentioning that these diagrammatic 
expansions are only obtained when the one-body, zeroth-order Hamiltonian is a diagonal 
operator, permitting the left side of Eq. pij) to be evaluated to give terms containing factors 
of orbital-energy differences, e.g., (e* — ^D- 

Note that Eq. (ji^ can also be expressed by 



q{n) TT^ 



V /en 



(45) 



where the reference state projector is 



(46) 



Expression ()45|1 is easily obtained from Eq. ()44|1 by using the definition of an open operator. 



adding (P$ 



P$ = 0), and using the following identity: 



q(n) rrq, 



P$ — —Q<^ (Hq S^^^\ P$. 

V / cn 



(47) 



D. Rayleigh-Schrodinger perturbation theory and the Bloch Equations 



Using Eqs. ©, (fTUj) . and pUj) . the Schrodinger Eq. ([T)) can be written 
The operator form of this Eq. is the Bloch Eq. j^, 74 . Iscl Isij] : 



i/r2ij>p$ — r2$p<j)i^f f2<j)P<j) , 



(48) 



(49) 



which yields Eq. (I48j) when multiplied by |$) from the right side. Substituting Eq. (p^ into 
and u sing and (jSE)), where i^o is Hermitian, gives a variant of the Bloch equation 



6^ 
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(£■0 - Ho) n^p^ = Q$ {vn^ - n^p^vn^) p$, 

where we have used the following identity: 



(50) 



p$ {vn^ - n^p^vn^) = o, 



(51) 



and this identity follows from intermediate normalization, given by Eq. Q, i.e., (P$n$P$ = 

For Rayleigh-Schrodinger perturbation theory, the wave-operator is given by the 
order-by-order expa,rs,o,r mB- Substituting this expression mto Eq, EB and equating 
the individual orders, gives |56|, [7^ 

(52) 

where the second term on the right side does not appear for (n = 1). 



n-l 
m=l 



E. The Linked diagram theorem 

When the zeroth-order Hamiltonian is in the diagonal, one-body form, give n by Eg. (jll]), 
it can be shown that the wave-operator il^ satisfies a linked diagram theorem 

(^0 - ^0*) ^1-^* = ^<i., (53a) 



where the individual orders, defined by Eq. ()33b|) . satisfy 

{Eo - H^) n^i^P^ = g$ (53b) 

and the additional / subscripts indicate that only the linked portions contribute - all dis- 
connected terms are open. 

In order to solve either of the above two Eqs, the wave operator Q<p is written as a sum 
of one-, two-, and higher-body excitations, 

Q^ = l + Qf + Ql + + ■■■ , (54) 

where 

1^-^ — ^ ^ ^y^j ^11? ; (55a) 



2 2! 



^ ^^wsx^lo-Wanj, (55b) 



rwsx 



^ 3! 



x^^g^fyalalalayaxCiw: (55c) 



rwsxty 



and Eqs. © remain valid. In addition, we require our coefficients to have exchange symme- 
try: 

■^rwsx ~ -^sxrun (56a) 
"^rwsxty '^sxrwty "^rwtysx^ .... (56b) 

The wave operator and its n-body operators fi* are also invariant to a unitary transfor- 
mation of either its occupied orbitals {ipo^ ^} or its virtual orbitals {ipu-^ 

The linked diagram theorem provides an exponential form for the wave operator Q^p] 
comparing the two forms of the wave operator, Eqs. ()24|) and (jSU), the following identities 
are easily proven: 

nf = Sf, (57a) 

nt = s! + ^{sf)\ (57b) 

fif = S! + SfSf + ^{Sf)\ (57c) 



IV. BRILLOUIN-BRUECKNER CONDITION 



Consider the Slater determinantal state, say |B), that satisfies the Brillouin-Brueckner 



condition 



26 



26 



38, 



for any single excitation from |0): 



(e;|i/|vi/) = o. 



|9^) = 4a^|9), 



(58) 



(59) 



where both the occupied and virtual orbitals determine the Brueckner determinantal- 
state |0): 



i^w e {ipo ©}, 



(60a) 
(60b) 



Using Eqs. ((T)) and (j^ . it is easily demonstrated that the wavefunction |\E') contains no 
single excitations from |0): 

1 



S 



(e;|i/|v[/) = (euxi/) = o. 



(61) 



Since the states |9^) are linearly independent, the wavefunction satisfies the following con- 
dition: 

F®|^) = 0. (62) 
where the projector for the singly-excited states is 

Pn= J2 E (63) 

and this subspace is completely determined by |0); is also invariant to a unitary trans- 
formation of occupied, or virtual, orbitals |3o|. 

Using Eqs. dHU), (jSS) and (jU^ . Eq. (j3H)) can be generalized: 

P^,H{1-P^,)\^) = 0. (64) 

The occupied set of orbitals {V^o^ ©} that satisfy Eq. (jMj) are called Brueckner orbitals. 
However, since these orbitals are invariant to a unitary transformation, Eq. ()64|) actually 
defines the Brueckner-determinantal state |G), since |0) determines P^®. 
Note also the following identities: 

nf = 0, (65a) 
Sf = 0. (65b) 

The first identity is obtained by substituting Eq. (0) into (jH^ and using and for 
($ = 9). The second identity uses either Eq. ()57a|) or Eqs. and (j^ . 
Since Sf is zero, we have 

fi© = e-^?ne. (66) 
Multiplying this equation from the right by |6) and using Eq. (0) gives 

|^) = e-^?|^). (67) 

Substituting this equation into Eq. (j58j) . and using Eq. (j63p . gives 

Pi^i7e-^?|^) = 0. (68) 
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This equation is the Brillouin-Brueckner condition for coupled cluster theory 
in Eq. (j64p . the Brueckner orbitals that satisfy Eq. (j68p are invariant to a unitary trans- 
formation, so Eq. (|UHj) defines the determinantal state |9), since |9) determines P^® and 
Sf. 



V. TRIAL WAVEFUNCTIONS AND ENERGY FUNCTIONALS 
A. General Requirements 

Consider four trial wavef unctions, denoted by where 77 = I, II, III, and IV. Each of 

these four states depends on the reference state |$), satisfies intermediate normahzation, 

($1^^^') = 1, (69) 

has no components within the singly-excited subspace, 

\¥i^) = {l-P,\)\¥i'), (70) 

and yields the exact state of interest when |$) is the Brueckner determinantal-state: 

l^g') = 1^)- (71) 
From these trial wavefunctions we can construct energy functionals: 

^,[<l>] = {^H\¥^') = + E^-^m, (72) 

where the correlation (co) energy-functionals are given by 

E^^l^m = ($ |//|vl/g^), (73) 

the trial correlation-functions are 

Q^in^) = Kl), (74) 

and Ei[^] is given by Eq. (fTT|. Operating on Eq. (f7T| by Qe and using Eqs. (fT3|) and 
we have 

I^Se) = I^Qe)- (75) 
Let us also define exchange-correlation (xc) energy-functionals: 

E^^m = E^^m-E^m, (76) 



where the exchange energy i?x['^'] is given by Eq. 

Eqs. (llUj) . (j71|) . and (j72p indicate that the energy functionals -E'r?['^'] yield the exact energy 
S when the reference state |$) is the Brueckner determinantal state |0): 

£ = E,[Q], (77) 



and from Eqs. (dH), (f7T|) . (fTSjl . and (fTUj) . the following identities are obtained for the 

correlation and exchange-correlation energies, £^co[0] and £^xc[0]: 

£eo[e] = E^^^[e], (78) 

^xc[e] = ^^,'[9]. (79) 

Substituting Eq. (HH) into Eq. (jEH), and using Eq. (HOI) for (|$) = |9)), gives the Brillouin- 
Brueckner condition for our trial wavefunctions: 

P®i7|^g') = 0. (80) 

We now define the explicit forms of these trial wavefunctions and give expressions for 
their correlation-energy functionals. Additional expressions for the correlation-functionals 
in terms of their wave-operator, or cluster-operator, amplitudes, e.g., tf^g^, are given in 
Sec. IVI| diagrammatic expressions are presented in Sec. IVIII (^1 



B. The first trial wavefunction 

The first trial-wavefunction is given by 

l^^') = (l-At)l*)- (81) 

It follows from Eqs. (jH^ and that Qf exclusively generates the singly-excited portion 
of the orthogonal space: 

p0fif|$) = ^]f|$), (82a) 
P®1]^|$) = 0, n^l. (82b) 

Therefore, this trial wavefunction can be written 

|vl/«) = (fi^-fif)|$). (83) 

Using this expression and Eq. (fT9|) . after substituting Eq. (f74j) into (f73|l . yields the first 
correlation-energy functional: 

Em=[H{X'^-nt)]^^, (84) 



where we also used (P$(x4> ~ ^f)!"^) = 0)- 



C. The second trial wavefunction 

The second trial-wavefunction is given by 

I^W) = e-^?|^). (85) 
Using Eqs. (0) and ()24a|) . this Eq. becomes 

|^m) = e(^*-^*)|<l>), (86) 
where we used the identity, given by 

e-^*n^ = e(^*-^?), (87) 



and this relation follows from Eqs. since Sq, and Sf commute j85|. 

Using Eq. (jHBj) after substituting Eq. (fflj) into (fTSj) . gives the second correlation-energy 
functional: 

where we also used the following: 

($|iJF$e(^*'^*)|$) - = 0. (89) 

D. The third trial wavefunction 

The third trial-wavefunction can be generated by its wave-operator: 

= |vl/D, (90) 

that can be expressed in an exponential form, 

Cl^ = e^^ = l + S^ + ^Sl + ^Sl + -- - , (91a) 

where 5"$ can be written as of sum n-body excitations, with the exclusion of a one-body 
operator: 

S^ = + Sf + ■■■ . (91b) 



The individual amplitudes are defined by the following equations: 



-y 

91 



t t 



21 / ^ ^rwsx'^r'^s^x'^w^ 



- y 



t t t 

rwsxty 



(92a) 
(92b) 



Let us also mention that the orbital convention, Eqs. ©, remains valid; we also require the 
coefficients to possess exchange symmetry, as in Eqs. (j^H))- 

Using Eqs. (PU)) and (j91ap . after substituting Eq. (fTij) into (f7^ . gives the third correlation- 
energy functional: 



H I e^* - 1 



where we used the following: 



($|i/P$e^*|$) - = 0. 

We define as a solution to the following variant of Eq. (j29|) : 



(1 - p*) (nn^ 



0, 



op,cn 



(93) 



(94) 



(95) 



which defines the trial functional I^E'^"') using Eq. (jHO)). 

As in Eq. (jHUj) . the operator HQ<^ can also be written as a sum of zero-, one-, two- and 
higher-body excitations: 



where we use the identity, given by 

and this relation follow from Eqs. (f?^ and (|^. (See also Appendix iBl) 
Substituting Eq. ^ into Eq. ^ yields 



(96) 



(97) 



HQ. 



op,cn 



HQ, 



0. 



(98) 



op,cn 



Substituting expansion into Eq. ()98|) and noting that each term is linearly independent, 
we have 



HQ, 



0, n>2, 



(99) 



op,cn 



and this relation can be used to obtain the coupled cluster equations for the amplitudes. 

We now demonstrate that |\&^"') is a valid trial wavefunction: Eq. (|7T| and the other 
relations for Sec. ()V A|) are satisfied. 

Consider a determinantal state, say |0'), that we require to satisfy the following condition: 



P®'(//(i©/) =0. (100a) 

V / op,cn 

Using Eq. (jHEI), we can easily verify that the following conditions causes Eq. ()l()()a|l to be 
satisfied: 



0. (100b) 



op,cn 



Adding Eqs. and (jlOOaj) . for (|$) = |6')), and comparing the result with Eq. 
indicates that 

ne' = ne,, (101) 

and therefore we have 

(Combining Eqs. (jUUj) and ()100b|) . and comparing the result with Eq. (jH^ . also yields the 
above two Eqs.) 

Multiplying Eq. (fTTl^ by Pf^' and using Eqs. and (jnT|) gives 

Pif|^) = 0. (103) 

Comparing this Eq. with Eq. ()(j2|l indicate that |G') is the Brueckner state: 

|0') = |0). (104) 

Substituting this result into Eq. ()102|1 indicates that Eq. (f7T|) is satisfies. All other relations 
from Sec. ()V A|l are easily verified. For example, Eqs. (|77j) . (f7H|) . and (jZHI) follow from 
Eqs. (IZ21) through (f7B|l . 

Substituting Eq. ()104|) into Eqs. pOO|) yields the following identities: 

Pi^ (Hne) = 0, (105a) 

V / op,cn 



= 0. (105b) 

op,cn 



These equivalent relations are alternative representations of the Brillouin-Brueckner condi- 
tion, since if they satisfied, then Eq. (jHUj) is also satisfied. Note that n$ does not possess a 
single-excitation operator, i.e., Sf is absent in Eq. ()91b|) . 

Using the definition of an open operator, Eqs. (j95p and (jl(J5ap can be represented by the 
following: 

g$ (Hn^) P$ = 0, (106) 

V / cn 

P® (H{le) Pe = 0, (107) 

\ / cn 

where 

Q$ = Q$-Pn. (108) 

As in Eqs. fj33|) . the wave- and cluster-operators of interest are given by order- by-order 
expansions: 

= + 4'^ + 4'^ + • • • > (109a) 

= fij) + f^J) + + . . . , (109b) 

where (fi^'' = 1.) Substituting these relations into Eq. ()91a|) and equating each order, we 
get similar identities as in Eqs. pljl : 

= 4^ (110a) 

1 

2 



= ss^+UsS'Y^ (iiob) 



and so so. 

From Eq. (jHSj), and using Eqs. (jl^ and p09|) - as in the derivation of Eq. ()44p - we have 
- (1 - a!) K4"0 = (1 - ^n) "'0 ■ (111) 

\ / op,cn \ / op,cn 

Solving Eqs. (jlll|) and (jllO|) in an iterative and sequential manner generates a linked diagram 
expansion for and a connected expansion for S^. The diagrams representing E^™^ [$] and 
4"'' are discussed in Sec. IVlIIBl and elsewhere j57|, respectively. 

Of the four trial wavefunctions I^E'^'), the third one |\E'^"'), we believe, is the most appli- 
cable; the fourth one is presented below for completeness. 



E. The fourth trial wavefunction 



The fourth trial wavefunction is a solution of the Schodinger Eq. within the subspace 
that neglects the single-excited states: 

(1 - FfO Hl^D = i?iv[$]|^r)- (112) 

From the variational theorem, it follows that the above energy functional provides an upper 
bound to the exact energy: 

Eiv[$] > S. (113) 

We now prove that the exact wavefunction and energy satisfy Eqs. (f7T|) and (f77|). where 
(r/ = IV). 

The proof uses the Schrodinger Eq. which can be written 

+ (1 - P®) H\^) = SP^^l^) +S{1~ P®) 1^), (114) 

where he have added and subtracted P® terms. Eqs. (j3Hj) . (jHH), and (jU^ indicate that the 
first terms on the right and left sides vanish, so we have 

(1-P®)^f|vl/)=^(1-P®)|^). (115) 

Eqs. (f7T|) and (f77|) . for (77 = IV), are obtained by comparing Eqs. (|112|) and (jll5|) . and using 
(j7(J|) . All other relations from Sec. ()V A|) are easily verified. 

As in the exact wavefunction of interest the trial wavefunction can be gener- 

ated by a wave operator f2$: 

n^\<!>) = |vl/g-'), (116) 

where is similar to - defined by Eqs. and - except that there is no excitation 
operator into P*: 

= 1 + nf + nf + ■ ■■ , (117) 

where 

^2 = ^ Yl ^twsxO'lo'Wan^, (118a) 

rwsx 

^3 = ^ XI ^fwsxtyaWsalaya^a^, (118b) 

rwsxty 



Using Eq. ()116p . after substituting Eq. (ffi|) into (fTHj) gives the fourth correlation-energy 
functional: 



cl 



(119) 



where we added (($|i7F$(^^$ - 1)|$) = 0). 

It is easily proven that is a solution to the following variants of Eqs. ()33b|l and ()52|1 : 



n$ = 1 + n^^^ + fii'^ + 



n-l 



m=l 



Pa 



(120) 
(121) 



where Q$ is defined by Eq. fll08|) . and we require Hq to satisfy 

P*Fo (1 - Ffi) = 0. 



;i22) 



VI. EXPRESSIONS FOR THE CORRELATION-ENERGY FUNCTIONALS 



Using Wick's theorem 
and two-body parts: 

where 



the Hamiltonian can be separated in to zero-, one-, 

(123) 



H = Hf + Hf + Hf, 



Hf 



J2^i\h\j]{alaj}, 



ijkl 



where the Fock operator is given by 



(124a) 
(124b) 

(124c) 



(125) 



and the identities within Eqs. p24a|l follow from Eqs. (jllll and (fTTj) . 

Substituting Eqs. (fT^ and p23|) into (fTHj). and using (jMj) . yield an expression for the 
correlation energy: 

£,M = {HM + H!nt)^^, (126) 



and from Eqs. ()124|) and ()55p we have 

^co[$] = ^^xfjwlF^lr] + IY1 ^rwsx {[wr\xs] - [ws\xr]) , (127) 



2 

rwsx 



where the coefficients are are assumed to have exchange symmetry, defined by Eq. ()56p . 

As is the correlation-energy expression (fTH|) . the correlation-energy functionals, given by 
Eq. (f7H|). can be written 

Eil^[^] = {Hxl),,, (128) 
where the trial correlation-operators generate the trial correlation-functions: 

Xim = (129) 

and I^q'j,) is given by Eq. (j74|) . Comparing this definition with Eqs. (j72j) and (j73|) . and 
using Eq. ()124a|) . we have 

E,[^] = (1 + xl) m = [Hil + xDL ■ (130) 

Comparing Eqs. (jSlD, (|HH1), (1221), and (fTT!I|) with (fnH|) . yields the following relations: 

xl = X* - , (131a) 

X'^' = e^* - 1, (131c) 
X'^ = n^-l. (131d) 

Substituting Eqs. (HH) and (fT^ into ((THI) for (|$) = |e)), we get 

xl = Xe, (132) 

and this expression indicates that any of the trial correlation-functions - Xe, Xe, Xe ' 
Xe ^ c^'^ t)e used to obtain the Brueckner one, xe- 

Substituting Eqs. (HSH) into Eq. (fnH|l . and using Eqs. (HD), (ED, (EQ), (ITT7I) . and 

give expressions for the correlation energy functionals: 

E»] = {Hfnt)^,, (133a) 
Km = {HfSf)^^, (133b) 
= (133c) 

^r[*] = {h^^i)^,- (133d) 



Using Eqs. ()55b|l . ()25b|) . ()92a|) . pi8ap . and ()124c|) . the above energy functional above can 
be evaluated, giving the following relations: 

rwsx 
rwsx 
rwsx 

E^Tm = ^ E ^rwsx i[wr\xs] - [ws\xr]) , (134d) 

rwsx 

where the coefficients are assumed to satisfy exchange symmetry, e.g., Eqs. (j26|l and (j56|) . 

For later use, we also mention that the correlation energy and the first correlation-energy 
functional, given by Eqs. (jl26|) and (jl33ap . can be written using the amplitudes: 

(135) 

E»] = (^iHfSfSf + HfS!^ ^, (136) 
where we have used Eqs. (|57ap and ()57b|) . 



VII. EXACT FOCK OPERATORS 

Consider generalized, or exact, Fock operators that are defined, in part, by the 

following matrix elements: 

(H^i^'li^w) = {<l>'jH\<il^'); r] = l, II, III, IV, (137) 

where the w and r orbitals are occupied and unoccupied within |$), respectively, as noted 
by Eqs. ©. By multiplying Eq. (jHUj) from the left by (0^|, using Eq. (jUH)) . and comparing 
the resulting relation to the above Eq., we have 

{iJr\:Fl''\^w) =0, (138) 

where the orbitals are defined by Eq. (j6(J|) . When satisfied by all orbitals, this expression 
is equivalent to the Brillouin-Brueckner condition, given by Eq. (|8(J|) . The operator form of 
Eq. (HHHl) is 

ke:F^''% = 0, (139) 



where 7* is the one-particle, density- matrix operator for the determinantal-state |$) j8| 



87 



88, 



89|: 



%= Yl (140) 

a:e{i/'o— »*} 

K$ is the projector into the virtual-orbital subspace: 

K$ = \i^r){A\, (141) 

and the identity operator / can be expressed by 

i = % + k^. (142) 

Multiplying Eq. (jl39p from the left and right by (-0^1 and {i/Jw) gives Eq. (jl38|) . 

Since all of our generalized Fock operators - Jq \ ^e^\ ^e^^\ and J^'^' - satisfy Eq. ()139p . 
any one can be used to define an exact Fock operator Pq: 

keTele = «^e-^e ^Te = «^e-^e ^Te, (143) 
and the Brillouin-Brueckner condition, Eq. p39p . becomes 

kePele = 0- (144) 
Using Eq. ()142|) . this Eq. can be written as 

(/ - 7e) -^676 = 0. (145) 

Since 7e is idempotent, 

7e7e = 7e, (146) 

Eq. ()145|) can be written as 



^^e7e - 7e^ej 7e = 0. (147) 
By requiring J-q to be, at least in part, Hermitian: 

%:Feke = 0, (148) 

yields the following identity: 

^e7e - 7e^e ) = 0. (149) 



Adding together Eqs. (jl47|) and (jl49|) . and using Eq. (jl42p . indicates that and J^q com- 
mute: 

^e,7e]=0. (150) 

Eq. (llSUp is a generahzation of the one obtained for Hartree-Fock theory |50|. 

Note that for any reference state, say |$'), we can find a corresponding state, |$), in 
which the following relation is satisfied: 

k^H"^^^ = 0. (151) 

Solving this expression in an iterative and self-consistent-field manner leads to the Brillouin- 
Brueckner condition, Eq. ()144|) . being satisfied, since when (|$) = 1'^'')), we have (1$) = |0)). 
Consider now the following application of the identity operator: 



where 



e 



Setting ($ = 6), and using Eq. p38p . gives exact Hartree-Fock Eqs: 

where the orbital lipw) is from the set {ipo — > ©}• 

Returning to Eq. p37|) . inserting the identity operator - defined by Eq. (fT^ - we have 

(AiH^'ii^^) = (Kim) + (KimQi), (155) 

where the trial-correlation functions |^q^,) are given by Eq. (f?^ : the first term on the right 
side of Eq. ()155|1 is the off-diagonal block of the Fock-operator F$: 



(Aihii^n.) = {'^i\Hm, (156) 

where is F$ defined by Eq. p25p . Substituting Eqs. p56|) and (jl29p into ()155p . we have 

(AlH^'li^^n) = (Alhli^n.) + {'^l\Hxm (157) 



where we use the more restrictive definition of an open (op) operator, presented in Ap- 
pendix |B| 

Setting ($ = 6), and using Eqs. (|132j) and (|138|) . gives another variant of the Brillouin- 
Brueckner condition: 

(Pe + iHxe),) =0, (158) 

V /op 

where this expression acts within the one-body sector of the Hilbert space, even though the 
subscript op indicates the open portion - defined by the n-body sector. 

Now let the second term on the right side of Eq. ()155|) define the off-diagonal block of 
correlation potentials v^^ix), given by 

{Alv'M..) = {K\H\¥^l). (159) 
Similarly, exchange-correlation potentials v^^{yi) are defined, in part, by 

(AKM = (^rlCIV^^-) - (AlK^l^u.). (160) 

Using Eqs. ()155p . ()156|) . and ()159|) . and with no loss of generality, our exact Fock operators 
J^''' can be written 

^^'" = F<, + C. (161) 

Multiplying Eq. ()161|) from the left and right by {iprl and {ipw), and using the one-body 
partitioning method of Eqs. ()B8|) . gives 

Comparing Eq. (jl57j) with (jl62|) . and using ()B9|) . we have 

(Oex= [(^X^)l]op' (163) 

and the Brillouin-Brueckner condition ()158p becomes 

f^e) = 0, (164) 

V /ex 

where 

Te = Fe + vZ, (165) 
and the 1] superscript is suppressed, since, in general we have 

(Oex = (Oex=(C)e.- (166) 



The remaining matrix elements of v"^^ - [i'w\v*o\i'x\ and (■^rl'^^co'lV's] ^ are at our disposal. 
By defining these matrix elements in a manner that is independent of i], but dependent on 
1$), f®o(x) and are completely, and unambiguously determined; our exact Fock operator 
can be diagonalized: 

^eV^f(x)=5f^f(x), (167) 



where orbital energies ef can be defined to 



give exact ionization potentials and electron 



affinities - exact Koopman's theorems 2l|, |57|. In addition, since the operators, /* and 
/*, that give the zeroth-order Hamiltonian Hq , and the exact Fock operators, !p^\ are not 
mutually exclusive, one tempting choice is 



t-i ■iri-n) 

Jo Ju $ ■ 



(168) 



By using the diagrammatic expansion for x^, and Eq. ()163|) . a diagrammatic expansion 
for (t'*o'')ex can be obtained that is a subset of the open one-body diagrams of [57] . 



VIII. REFERENCE STATE ONE PARTICLE DENSITY MATRIX THEORY 



A. Functionals of the one-particle density-matrix 7 



There is a one-to-one correspondence between the set of determinant states, {|$)}, and 
;heir o ne-p article density-matrices 3, Is^, {7}, where these density- matrices are given by 
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88, 



7(x,xO= ^-WC(x')- (169) 

Because of this correspondence, determinantal states are uniquely determined by their one- 
particle density-matrix: [$(7)); functionals, or functions, that depend on |$), can be written 
as ones depending on 7. For example, the total energy £^ Eq. (fTO|l . and our energy functionals 
Eq (IZ2I), can be written 



£ = £;i[7]+^co[7], 
E,Yi] = i?i[7] + i?S'[7], 



(170) 
(171) 



where, in addition, our trial wavefunctions |^^(^)) can be denoted by \^^^^). 



For simplicity, we require the external potential v{r) to be a spin-free operator, so the 
first-order energy can be written as 

Eib] = j [-iVr'7(x,x')]^,^^dx + j ^(r)7(x,x)dx + i?j[7] -i?x[7], (172) 
where the Coulomb and exchange energies are 

EAl] ^r2S(xi,Xi)7(x2,X2) C?XidX2, (173) 

EAl] ?'r2S(xi,X2)7(x2,xi) dXi(iX2. (174) 

(As in Eq. (0}, an integration over Xj actually implies a summation over the spin variable 
uji and an integration over the spatial portion r^.) 

Similarly, the one-body operators Eq. ()16ip . can be written 



P^^ = F, + vl2, (175) 
where, instead of Eq. ()125j) . the Fock operator is given by 

F^ = -lV'^ + v + J^-K^, (176) 
and the Coulomb J-y and exchange operators satisfy: 

J^0(xi) = j rf2S(x2,X2) 0(Xi) rfX2, (177) 

K^0(xi) = y"rr2W(xi,X2)0(x2)dx2. (178) 
In addition, the identity operator, given by Eq. ()142|) . can be written 



i = l + k^, (179) 
where the density-matrix operator 7 is defined by its kernel, 7(x, x') 8^: 

70(x) = j 7(x,x')0(x')rfx', (180) 
and fi;^(x, x') is the kernel of the virtual-space projector: 

/«^0(x) = y K^(x,x')0(x')(ix', (181) 



where 

«:^(x,x')= Yl M^Wri^'), (182) 

and kq, is given by Eq. p41|) . Similarly, the one-particle density-matrix for the Brueck- 
ner state, say r(x, x'), and its density-matrix operator, say f, are given by the following 
expressions: 

r(x,x') = Yl ^-(x)^l(x'), (183) 



f(/)(x) = y r(x,x')0(x')cix', (184) 

where we have 

|e) = |$(r)). (185) 
Since the one-particle density-matrix, f, also satisfies 

f = 7e, (186) 

where 7e is given by Eq. ()140j) . the Brillouin-Brueckner condition, given by Eq. ()144|) . and 
its complex conjugate, given by Eq. ()148|) . become 



krJ^'rT = 0, (187a) 
f:Frkr = 0; (187b) 



the commutation condition, given by Eq. ()15U|) . can be written 

^r,r] =0; 

the exact Hartree-Fock Eq. ()154|) is 



(188) 



x€{ipo-*T} 

where the occupied orbital, |^^), is from {ipo — > t}; furthermore, the other Brillouin- 
Brueckner condition, Eq. (jl64j) . can be written as 

(:Fr) = 0, (190) 

V /ex 

where Eq. ()165|) becomes 

:Fr = F^ + (191) 



In the following subsections, we illustrate how the correlation energy and correlation- 
energy functionals, £co[l\ and £"^^'[7], can be obtained from perturbation theory, in which 
all terms (or diagrams) explicitly depend on 7. (An explicit expression for the first-order 
energy, i?i[7], is given by Eq. (jl7!2j) .) 



B. The correlation energy as a functional of the one-particle density matrix: £co\l\ 

Consider the normal-ordered form of the Hamiltonian if, given by Eqs. ()12Hj] and ()124j] . 
where we choose to represent this operator using the following orbital sets: {tfjo'^ lif2} and 
{V'm^ iJz}^ and both sets are introduced in Sec. IIII CI Since the one-body portion of the 
Hamiltonian is determined by the Fock operator F^, it seems appropriate to use the 
following notation: 

Hi = {F^}, (192a) 

or, equivalently, the expression {F-y} is defined by the following procedure: Obtain the Fock 
operator F^ for the one-particle Hilbert space, given by Eq. (jl76j) : write this operator in its 
second quantized form using the true vacuum state | ); re- write using normal ordering 
with respect to the new vacuum state, |$(7)): (F^ = (-Fy)ci + {-^7}); the uncontracted term 
is IF, I . In other words, \F^} is the uncontracted term when the operator F, is written in 

nnnn 

normal ordered form using Wick's theorem |55l . 1561 16 IL I73| . 

Using this notation, we also write 

H^2 = {ri2h, (192b) 

where, again, is written in second quantization using the true vacuum state | ) - given 
by the second term on the right side of Eq. ^ - and then it is re-written using normal- 
ordering with respect to the new vacuum state; the uncontracted term is {rfg^}^, where the 
additional 7-subscript from {. . .}^, serves to remind us that the vacuum state is [$(7)). 
(This 7 subscript is suppressed in Eq. ()192a|) : The vacuum state is understood, since 7 
determines Fy.) 

The correlation energy Sco is determined by the one- and two-body parts of the Hamilto- 
nian, {Fy} and {rjfg^j^,. However, the individual orders of the perturbation expansion for £co, 
also depends, in addition, on the zeroth-order Hamiltonian, given by Eq. (j^T|) : this operator 



can be written as 

H^ = d^ + %, (193) 

where these terms - and il^ - are the occupied and unoccupied portions of Hq - (-f^o )oc 
and (-f^o)un ^ and are given by the following: 

il^ = ^ e^^alar. (194b) 

Using this notation, our zeroth-order Hamiltonian in normal-ordered form can be written 

as 

H^ = Eo[^] + {d^} + %, (195) 

where M-y is already normal-ordered; the constant term Eo['~f] is the zeroth-order energy of 
|$(7)): 

//J|*(7)) = £^o[7]|$(7)), (196) 

and is given by 

^o[7] = E (197) 

Note that the first-order and the correlation energies, Ei['~f] and Sco, do not depend the 
zeroth-order energy -Eo[7]- 

The perturbation V^, defined by Eqs. (I42j) . can also be written in normal-ordered form: 

V, = + + V^, (198) 

where, from Eqs. (|123|) . (jl24ap . (jl92p . and (|195|) . the individual terms are given by the 
following expressions: 

= i?i[7]-i?o[7], (199a) 

V7 = {F^} - {0^} - (199b) 

V7 = {r^,%. (199c) 

The one- and two-body parts of H are given by Eqs. ()124b|) and ()124c|) . and are equal to 
{F-y} and {r^^g^}^, as indicated by Eqs. ()192p . The Goldstone diagrammatic representation 



of these operators can be written in the following manner |y , |5a Isa ly, |5a |56|, IH ii 



m = m 



{'"l2^}7 



(200a) 



(200b) 



The one-body part of the perturbation Vi is usually represented by a single diagrammatic 
operator. However, for our purposes, it is convenient to use separate diagrammatic operators 
for the three terms on the right side of Eq. ()199b|) . where {F^} is presented by Eq. ()200ap . 
Since the other two terms are diagonal, it is appropriate is simply represent them as (unfilled) 
arrows: 



V 



A 



(201a) 



(201b) 



In contrast, hole- and particle-lines, by themselves, are represented by filled arrows: Vcind A 
As a slight alternative to the usual approach to evaluate the diagrams of the correlation 
energy Sco 3 Q, H Q H H H Q, 

we associate a hole line corresponding to a 
ty-occupied orbital with a '0^(xi)?/^^(x2) factor; we associate a particle line corresponding 
to an r-unoccupied orbital with a '?/'r(x2)'?/'^(^i) factor, where Xi and X2 denote the dummy 
integration variables that arise from the vertices. Using this convention, the sole diagram 
involving the Fock operator from second-order perturbation theory can be evaluated in 
the following manner: 




i^rw) ^ j rfXidX2 (^F^i7^(Xi,X2)j ■ F^27r(x2,Xi), 



where 



(202) 



(203) 



and the repeated indices - r and w - are summed over; F^i denotes the Fock operator F^ - 
given by Eq. ()176p - acting upon (xj); the term (F^^ ■ ■ ■ )■ indicates that F^i exclusively acts 



within the brackets; furthermore, the wth component of the (one-particle) density-matrix 7 
is denoted by 

7^(xi,X2) = V^^(xi)V'* (X2); (204a) 
the rth orthogonal-component of 7 is denoted by 

7,(xi,X2) = ^,(xi)^;(x2), (204b) 

where, for a complete set of orbital states, we have 

5(xi - X2) = ^7^(xi,X2) + ^7^(xi,X2), (205) 

w r 

which is a shorthand notations for 

5(xi - X2) = (5(ri - r2)5^ia.2- (206) 

In order to further compress our notation, we use the convention that all repeated dummy 
indices are integrated over and restrict the Fock operator F^^ to exclusively act upon the 
first index of any two-body function, i.e., (F^ia'(xj, Xj)a(xj, xj) = a'(xj, Xj)F^ia(xj, x^)); 
Eq. ()202j) can then be written as 



= i^rw) ^^7i7«.(xi,X2)F^27,.(x2,xi), (207a) 
and the other two diagrams from second-order perturbation theory have the following forms: 



^i^rwsx) Vi2V3/7^(xi,X3)7^(x3,Xi)7^(x2,X4)7.(x4,X2), (207b) 



-\{^rwsx) ^^12^^^3/7^«(Xl,X3)7^(x3,X2)7^(x2,X4)7,(x4,Xl), (207c) 



where 



'rwsx ^rw ~^ ^sx' (208) 



The diagonal terms arising from the zeroth-order Hamiltonian, given by —{67} and — W^, 
and represented by Eqs. ()201|1 . first appear in third order. For example, the following two 



diagrams can be obtained by inserting —{O^} and —il^ into the diagram on the left side of 
Eq. (jWifl : 



(209a) 



r^7l7t« (Xl , X2) F^27r (X2 , Xi ) 



(209b) 

The hole-hne operator {67} generates an additional hole line when inserted into a diagram 
and, therefore, a factor of —1 is included when diagram ()209a|l is evaluated, where this 
factor cancels the —1 factor from — e^. Since this type of cancellation always occurs, as an 
alternative, we associate a factor of for {O^} insertions, and treat vertices as ones 
that do not generate additional hole lines; is associated with a — factor. Keep in mind, 
also, that these operators generate an additional energy-denominator factor, e.g., ejl^, when 
inserted into a diagram. 

One advantage of partitioning the one-body part of the perturbation into individual 
components, as indicated by Eq. ()199b|) . is that it yields correlation-energy diagrams that 
explicitly depend on the Fock operator F^. Furthermore, an overall dependence of the 
correlation energy E^o on the one-particle density matrix 7 becomes, to a certain extent, 
transparent, by using Eqs. ()204|) . p69|) . and p82|) . yielding the following identities: 



7(xi,X2) = ^7^(xi,X2), (210) 
fi;^(xi,X2) = ^7r(xi,X2), (211) 

r 

and note that depends, exphcitly, on 7: 

5(xi - X2) = 7(xi,X2) + «;^(xi,X2), (212) 

where the operator form of this relation is Eq. ()179j) . 

The individual diagrams depend, in part, on each of the 7^ components, given by 
Eq. ()204ap . and the orthogonal components 7^, given by Eq. ()204b|) . In addition, each 
diagram depends on the set of orbital energies {e^}, which are at our disposal. In order to 



make each diagram an explicit functional of the one-particle density matrix 7, we choose all 
occupied orbitals to be degenerate, with energy el\ also, we choose all unoccupied orbitals 
to be degenerate, with energy e^. With these choices, the zeroth-order Hamiltonian, given 
by Eq. (jUJ, becomes 

HI = E + ^« E "J^- (213) 

and since this operator is invariant to a unitary transformation of occupied or unoccupied 
orbitals, it no longer depends on /J and - any set of orbitals defining 7 is appropriate - 
so we can write 

iil = ^lY. «t«- + e2E (214) 

It is easily proven that all perturbative orders, except for the zeroth-order, depend on the 
orbital-energy difference e"^, given by 

e, = el- el, (215) 

and not on the individual orbital-energies, el and e^. Therefore, we can choose (e^ = 0), 
and so our only parameter is e^. With this choice we have 

H^o = e^N^, (216) 

where iV^ is the number operator for the occupied orbitals, 

^7= E (21'^) 

■u)G{'0o^7} 

and it gives the total number of occupied orbitals when acting on a single determinant. In 
the one-particle Hilbert space, this operator is the projector for the occupied subspace - 
spanned by {ipo^ 7} ^ or, the one-particle density-matrix operator: 

^7= E l^-)(^-l=7- (218) 

Using the above two expressions, let us generalize the definition of 7: 

7= E (2^9) 



and write the zeroth-order Hamiltonian in a simplified form, given by 

i/J = ^7 7- (220) 

By normal-ordering this expression, we have 

H2 = e,N, + e,{^], (221) 

where A^-^ is the number of particles within |$(7)), and from Eq. ()195p . we get the following 
identities: 

EoU = e^N^, (222) 
{O,} = £,{7}, (223) 
= 0; (224) 

furthermore, our zero- and one-body portion of the perturbation, Eqs. ()199a|) and ()199b|) . 
have the following modified forms: 

1/^ = E,[^]-e,N„ (225a) 

^ (225b) 

Eq. ()224p indicates that the unoccupied operator, u^, represented by Eq. ()201b|) . does 
not appear in the expansion of the correlation-energy 8^0', {O^}, represented by Eq. ()201a|) 
and given by e-y{j}, is associated with a factor of Sy. Each diagram now becomes an explicit 
functional of 7 and k^. For example, the second-order diagrams can be written in the 
following manner: 



£:^F^i7(xi,X2)F^2«;7(x2,Xi), 



(226) 



^'^12^^34S(Xl, X3)/«^(X3, Xi)7(x2, X4)«;^(X4, X2), 



(227) 



4 



£7 ^'^12^?^3/7(X1, X3)/€.^(X3, X2)7(X2, X4)k^(x4, Xi), 



(228) 



where is given by Eq. 

In order to remove the exphcit dependence on k^, we first note that the kernel of the 

□ 

Fock operator is given by 

F^(xi,X2) =(5(x2-xi) [-iv,2 + v(r2)]+5(x2-xi) j r2-3S(x3,X3)rfx3-rf2S(xi,X2), (229) 
or we can use its matrix representation from a complete basis set: 

F^(xi,X2) = 5^^,(xi)(^,|F^|^,)V';(x2). (230) 

In any case, we use the following identity: 

F^ia(xi,X2) = y rfx3F^(xi,X3)a(x3,X2), (231) 



to modify Eq. ()226|) . yielding 

A jr = £:-1F^(xi,X3)7(x3,X2)F^(x2,X4)k^(x4,Xi). 
• V (232) 

Using Eq. ()212j) . this term becomes 

£::^^F^(xi,X3)7(x3,X2)F^(x2,X4)/t^(x4,Xi) = (233) 
e:^^ [F^(Xi, X3)7(X3, X2)F^(X2, X4)5(X4 - Xi) - F^(Xi, X3)7(X3, X2)F^(X2, X4)7(X4, Xi)] . 

The first term on the right side reduces to 

£:~^F^(Xi, X3)7(X3, X2)F^(X2, X4)5(X4 - Xi) = £:^^F^(Xi, X3)7(X3, X2)F^(X2, Xi), (234) 

where Fy(x2,xi) is given by Eq. ()23()j] . and not by Eq. ()229j] . since this form of the Fock 
kernel contains the laplacian V^, which is only defined when acting upon a function. As 
a possible alternative, the term on the left side of Eq ()234j) can be treated using the Fock 
kernel, given by Eq. ()229p . and a representation of delta function that is convenient to 
differentiate. This approach is probably more efficient, since it avoids the sums over the 
one-particle basis that appears in the matrix representation of the Fock kernel, Eq. (j23(Jj) . 

The substitution for k^(x4,xi) in Eq. (j233|) is applicable for any diagram. For example, 
consider the second order (correction to the) energy, denoted by -E'2[7], that is given by a 
sum of the right sides of Eqs. p32|l . p27jl . and (j228|l : where, using Eq. (j2T2|l . it can be 



E2[l] = e Ky(xi,X2)7(x2,X3)F^(x3,Xi) - F^(X1,X2)7(X2,X3)F^(X3,X4)7(X4,X1) 



(235) 



4^7 



ri2^7(xi,xi)7(x2,X2) - 2ri2Vi3V(xi,Xi)|7(x2,X3)|^ +ri2V3/|7(xi,X3)|2|7(x2,X4) 



- ^12 l7(xi,X2)| + 2r3i 7(^1, X2)7(x2, X3)7(x3, Xi) - r^^ r^^ 7(xi, X2)7(x2, X3)7(x3, X4)7(x4, xi) 

The individual terms from Eq. ()235|1 . as well as higher-order terms, can be represented 
by diagrams, where we write the Fock kernel in the following manner: 



and the two terms from Eq. 




F^(xi,X2) 



are given by 



-5 ^F^(xi, X2)7(X2, X3)F^(X3, X4)7(X4, Xi), 



(236) 



(237a) 




^F^(Xi, X2)7(X2, X3)F^(X3, X4)5(X4, Xi). 



(237b) 

These diagrams are evaluated by following the arrows in a backwards direction, where the 
circled arrow in Eq. ()237aj) indicates that this line is not a particle line, but a hole one with 
an additional factor of —1. This convention conforms to the one used for folded diagrams 
in valence-universal multireference perturbation theory 

0, y, Q, H Q- The 5 hue in 
Eq. ()237bj) yields a 5(x4,xi) factor; the integration over the X4 can be performed, as in 
Eq. ()234j) . and the resulting diagram can be represented by 




e F^(xi,X2)7(x2,X3)F^(x3,xi) 



(238) 



where the sole purpose of the above internal line without an arrow is to preserve the dummy 
index - in this case xi - for the two vertices it connects. 

The diagrams from Eqs. ()238p and ()237a|) correspond to the first two terms on the right 
side of Eq. ()235|) : the other terms are given by the diagrams from Fig. ^ The mirror-image 



FIG. 1: Second-order correlation-energy £co diagrams. 

diagrams (b) and (c) are equivalent, or non-distinct. Therefore, one of them can be omitted 
if a factor of 2 is included when evaluating the other one. Similarly, diagrams (f) and (g) 
are also non-distinct, so one can be omitted. Higher-order, correlation-energy E^o diagrams 
are obtained in a similar way. 

It is well known that £co is given by the set of connected diagrams^when fat least some 
of) the exclusion principle violating (EPV) diagrams are included 

0, 0|- Diagrams (a) through (f) in Fig. El present examples of diagrams that contribute to 
£co- (Diagram (g) is disconnected and does not contribute.) When orbital degeneracy is not 
imposed, the hole and particle lines in these diagrams correspond to 7^ and 7^, respectively. 
However, by imposing orbital degeneracy, as indicated by Eq. ()220|) . and using Eq. ()212p . 
these diagrams can be converted into ones that explicitly depend on 7, where the arrows 
representing the particle lines are either deleted or circled, in all unique ways. Additional 
diagrams containing the (or —{67}) interactions are generated by adding unfilled 

arrows, as represented by Eqs. ()201a|) and ()223|) . and including a factor of for each one 
appearing in a diagram, as well as the appropriate energy-denominators. For example, the 
following sixth-order diagram is generated from diagram (a) from Fig. |21 



£^^7(xi, X2)7(X2, X3)7(X3, X4)7(x4, X5)F^(x5, X6)7(x6, Xi), (^239) 

where the dummy index from the left-bottom vertex is Xi. The above diagram contains 
four hole-lines, one loop, and a circled particle-line, yielding, overall, a factor of +1. In 
addition, each one of the two unfilled-arrow interactions give a factor of e^,, and an energy- 
denominator factor of 2e^ - since both interactions appears within a double excitation - 
so overall a factor of {2e^)^^ obtained for each of these interactions; in general, 





FIG. 2: Examples of £co and Eco diagrams. 

an unfilled- arrow interaction gives a factor of ^ when appearing on a n-body excitation. 
In the above diagram, an additional factor of also appears from the remaining three 
energy-denominators. 

A discussion of perturbative convergence when imposing orbital degeneracies is presented 
in Appendix O 



C. The correlation-energy functionals -Eco [7] 

We now briefly discuss the diagrammatic representation of the correlation-energy func- 
tionals the details are given elsewhere 0|. The diagrams for E'^^J are a subset of the 
correlation-energy E^o diagrams. Therefore, we continue to use the diagrams from Fig. |21 
as examples, where, as mentioned previously, these diagrams can be converted into ones 
that explicitly depend on 7; additional diagrams containing —s^{^} interactions are easily 
generated from the ones in the figure. 

Comparing Eqs. and (jlHfij) . we see that the £co diagrams that contribute to {HlSl)c\i 
or ({F^}S'7)ci, do not contribute to i?™, where these diagrams have a final interaction, 
associated with if^, and a single-excitation below this last interaction, associated with SJ. 
Diagram (a) from Fig. |21is an example of a {HiSj)ci diagram. 

Eqs. (Unni) and (p3E|l indicate that the E^^^^ dia grams are a subset of the E^^ diagrams, 
where diagrams contributing to (|if2 "S*! 'S'7)ci do not contribute to E^^K Diagram (b) in 
Fig. 121 is an example of a {^H^S^SDci diagram when the top H2 interaction is 

removed, the resultant diagram is disconnected - it possessing two pieces, where, using 
the factorization theorem, it is easily demonstrated that each of these fragments contribute 
to Sj. In general, any i?^"' diagrams diagram that generates a disconnected diagram by 
removing its top interaction is not a i?^"' diagram, and the resulting disconnected diagram 
contributes to ^SjSj. 



It is easily shown that the i?co"' diagrams are a subset of the i?^"' diagrams [57[. In general, 
any i?^"' diagram that has an intermediate single-excitation does not contribute to i?co"'- 
An example is given by diagram (c) in the figure. This diagram possesses a single-excitation 
that appears in the second order. In other words, if we remove the top-two interactions, a 
second-order, single-excitation diagram is produced. 

Another type of i?^"' diagram that does not contribute to E^™^ is diagram (d). If we 
remove the top two interactions from this diagram we obtain a disconnected diagram com- 
prised of two fragments - one fragment being a third-order single- excitation diagram and the 
other one being a first-order double-excitation diagram. In general, if we remove any of the 
top interaction of a i?^"' diagram and obtain a disconnected diagram, in which one or more 

fragments are single-excitation diagrams, the (parent) i?^"^ diagram does not contribute to 
^(iii) 

CO 

Diagrams (e) and (f) contribute to E^^K Diagram (e) is connected in each order; diagram 
(f) is disconnected in the second order, but both fragments are double-excitation diagrams. 

By including certain EPV diagrams, disconnected diagram can be excluded from the set 
of diagrams that contribute to the first three correlation-energy functionals: E^^^, i?co', and 
Ej,™K However, disconnected diagrams, like diagram (g), do contribute to E^^^J^; the set of 
diagrams representing E^^J^ excludes diagrams with intermediate single-excitations, and this 
omission removes certain disconnected diagrams that are needed to invoke the factorization 
theorem and obtain a linked diagram expansion. 



IX. TREATMENT OF THE EXTERNAL POTENTIAL 
A. External potential dependence on £co and E'co 

Density functional theory employs a universal exchange-correlation functional, indepen- 
dent of the external potential, permitting approximations to be derived from model systems, 
where, in the vicinity of the model systems, the general form of the exchange-correlation 
functional is known. In contrasts, our correlation-energy functionals depend on the external 
potential v{r): £'^^'[7,1;], and are, therefore, in this sense, not universal. Nevertheless, by 
partitioning E^l^['j,v] into individual terms that are - to a varying degree - universal, we 
can also obtain approximations from model systems. 



In order to pursue this approach, and partition E^l^[y,v], we first partition the Fock- 
operator, given by Eq. ()176p . into two terms: 

= FF^ + V, (240) 

where F^ is the Fock operator for an electron gas: 

F^ = -IV^ + - K^, (241) 

and this operator is independent of the external potential v. Using this partitioning, the 
one-body part of the perturbation Vi , given by Eq. (fT99b|l . can be written as 

= {F^} + {v} - {O,} - U„ (242) 
where the normal-ordered, uncontracted portions of F^ and v have the following forms: 

{F^} = T.l'\^y\M4a,}, (243) 

ij 

{v} = J2[t\v\j]{ala,}. (244) 

The correlation energy and our correlation-energy functionals, £^co[7, v] and E^l^['y, v], have 
a dependence on the external potential v that arises, exclusively, from the {v} contribution 
to Vi , as indicated in Eq. ()242|) . In the diagrammatic treatment presented in Sec. IVIII B| 

is partitioned into three terms, given by Eq. (ITM . We now consider an approach where 
we partitioned it into the four terms, given by Eq. ()242|1 . 

The diagrammatic representations of —{O-y} and —{il-y} are given by Eqs. (jSHU); the 
other two terms from Eq. ()242|1 use the following: 



{m = o 



M = o 



(245a) 



(245b) 

For simplicity, in this section, we consider diagrammatic examples that do not distinguish 
between the various functionals, do not contain —{O-y} and —{il^} insert, and do not explic- 
itly depend on 7. For example, diagrams (a) through (d) in Fig. 01 are simple examples of 




(a) (b) (c) (d) (e) 



FIG. 3: Fourth-order E^o and Eco diagrams. 

fourth-order diagrams that contain the {F^} and {v} operators. These diagrams contribute 
to the correlation energy £co[l,v\, and any of the correlation-energy functionals E''^^[^,v\. 
Diagram (e) is the corresponding diagram that uses the {F^} operator, given by Eq. (j2(J(Jaj) . 
All four diagram, (a) through (d), can be obtained from diagram (e) by replacing the {F^} 
operators with the {F-^} and {v} operators in all unique ways; the sum of diagrams (a) 
through (d) is equal to diagram (e). In general, all diagrams containing {F-^} and {v} can 
be obtained from the {F^} diagrams by usin g th is approach. A similar procedure can be 
used when adding an additional perturbation j56l |. 



B. Electron gas terms 

Note that the diagrams from Fig. |S1 can be partitioned into one that does not depend on 
the external potential diagram (a) in the figure, and the remaining ones that do depend on 
V. This is a general result when using the partitioned Fock operator F^, given by Eq. ()24()j) : 
therefore, the correlation energy £co[^,v\, and the correlation-energy functionals E'^^[^,v\, 
can be divided in the following manner: 

^cobM = ^:r[7]+^^r'[7,^], (246a) 
E^Z' [l,v] = i?<r' [7] + E^c^^' [7, ^] , (246b) 

where i^co [7] -^co°'[7] universal functionals of 7, independent of the external po- 
tential; they are given by the sum of diagrams that do not contain {v} inserts; including 
diagram (a) in the figure, and the second-order Coulomb and exchange diagrams, given by 
Eqs. fl207b|) and ()207c|) . The remaining diagrams of Sco['l,v] and E^l^['j,v] contribute to 
£^^o+'[7,f] and E^l-^^''['j,v], respectively, including diagrams (b) through (d) in the figure. 
We refer to ^^co W -^co^'W electron-gas terms, since, for an electron gas, v is 



constant, so {v} is zero, and we have the following: 

^coW = 1 ^ ^ (247a) 

> (electron gas) 

E^l'il] = i?r'[7].J (247b) 

From Eq. ()78|1 . the terms on the left sides of Eqs. ()247j] are equal when 7 is the one-particle 
density-matrix for the Brueckner state; hence, we have 

^rN=^ry, (248) 

where Tg is the Brueckner-state one-particle density-matrix for an electron gas. For future 
use, we denote the correlation energy of an electron gas by £^^^""-'[7]; explicitly, we have 

£^^as) ^ y| ^ y ^ constant , (249) 
so, using Eq. ()247a|) . Eq. ()248p can be written as 

^rN = ^r'N. (250) 

C. Atomic and molecular terms 

In order to further generalize Eqs. ()246|) . we partition the external potential into its 
individual components: 

v = J2''m, (251) 

m 

where Vmi'r) is the external potential from the mth nuclear point-charge at R^: 

= - .^^"l (252) 
and we choose the following diagrammatic representation for this operator: 

{Vm} = 

(253) 

Using the decomposition of v, given by Eq. 1)2511) . the one-body part of the perturbation 
V^, Eq. (I212I), becomes 

= {An + E{^-} - iO,} - {U,}, (254) 




(a) (b) (c) (d) (e) 



FIG. 4: 8co and Eco diagrams with v components, 
permitting us to generalize Eqs. (|246p in the following manner: 

^.obM = + E^co'[7,^J +^:r[7,^], (255a) 

m 

Eil' b,v]= E^^^ [7] + ^rn] + E^^-^ (255b) 

m 

where the atomic terms, £^^^'[7, v^] and i?^2'^'[7) "^m], depend on they are the sum of all 
diagrams that contain one or more {vm} inserts, but have no inserts from other external- 
potential components, say {fn}, where {n 7^ m). For example, diagrams (b), (c), and (d) 
from Fig. m contribute to £co[l^ ^rn\ and E'^^^^b^ Vm\- These diagrams are generated from the 
corresponding diagrams in Fig. El by replacing all the external-potential inserts, {w}, given 
by Eq. ()245b|) . with the individual components, {fm}, given by Eq. ()253p . 

The remaining diagrams contribute to E^^^'b^v] and i?^^'^+'[7, v]. These diagrams have 
two or more inserts from (iz/ferent external-potential components, e.g., diagrams with both 
{vm} and {vn} inserts, where (m 7^ n). Diagram (a) in the figure is an example of a £^co^''[7, v] 
or E'^^'^^'ITi ^] diagram. Diagram (e) is also a i?^o^+'[7, v\ diagram, unless {m = n = p), where 
it then contributes to Eco['J,Vm] and £'^^'^'[7, f^]- 

For an atomic system we only have a single nucleus, and we get 

^cob,^^] = £M+£^b,Vn^, 1 (256a) 

> (atomic systems) 
Eil^b, v] = [7] + i?^r'[7, Vrn]. J (256b) 

The atomic terms, S^l^^bjVm] and E^l''^^['j,Vm], are - for a range of systems - universal 
functionals: They contribute to Sco[7,v] and E^l^b,v] for any system with a nuclei with 
a charge of Z^- Note that two atomic tems, say S^'^^['j,Vm] and i^^^' [7, are the same 
functionals when they are from the same nuclei {Zm = Zm'), but they are expressed from 
different nuclear locations, since we have (R^ 7^ R-m')- 



Factoring Zm from the diagrams, we get the following expansions: 

C [7, v.^\ = Z^S^l^l^ W + Z^S^oU W + ^^^3) W + • • • , (257a) 
' [7, vr.] = Z„,E^^^ [7] + Z^E^-'l^^ [7] + ZiE^lf^^ [7] + . . . , (257b) 

where i^cowb] and E'co'wM are universal functionals, independent of Zm] they are given 
by the sum of all diagrams with i inserts of {vm}, with {Zm = !)• For example, diagrams 
diagrams (b) and (c) from Fig. HI contribute to i^cowH and i?co(i)[7]) for {Zm = 1), since they 
have a single insert; diagram (d) contributes to i^co(2)[7] and £'co{2)[7]j diagram (e) contributes 
to £co(3)[l] and £^co(3)[7], ii {n = p = m). 

Further generalizing Eqs. ()255|1 . we have the following external-potential expansions: 

S,o[^,v] = £:r[7] + J2^cJ[l,Vm] + J2 ^cl'b,vm,vn] + • • • , (258a) 

m m>n 

Ki'hv] = [7] + E^cr'[7,^j + E Eir'b,vm,v^] + ..., (258b) 

m m>n 

where the diatomic terms, £^^^'[7, v^, and E'^l''^^['j,Vm,Vn], depend on Vm and Vn', they are 
the sum of all diagrams that contain one or more {vm} and {vn} inserts, but have no inserts 
from other external-potential components, say {vp}, where {p 7^ m) and {p ^ n). They are 
universal functionals for any system with one or more nuclei with charges of Zm and Zn, 
but, in addition, they also depend on the molecular geometry - the distance between the 
charges. Diagrams (a) from Fig. 0] is one example. By factoring the Zm and Z„ terms from 
the diagrams, analogous expansions, as in Eqs. ()257p . can also be obtained. Triatomic and 
molecular terms are defined in a similar manner. 



X. APPROXIMATIONS 

The diagrammatic methods presented in Sec. I Villi and in the previous section, gives 
explicit expansions for the functionals £co{li 'v] and i?co'[7> and for the terms on the right 
side of Eq. (j258|) . However, in the approximations we consider below, we often assume 
that these functionals can be expressed in a more simplified form, for example, as integrals 
involving the coordinates of only two electrons: 

^^cob,''^] ^co(xi,X2) (ixi<ix2, (259a) 

^^0^7,^^] = \jj G'i^o'(xi,X2)rfxidx2, (259b) 



where the integrands, ^co(xi,X2) and (^^^^(xi, X2), exphcitly depends upon 7(xi,X2), 
7(x2,xi), v{ri), and v{r2), and can include gradients or higher-order derivatives, e.g., 

Vh(Xl,X2). 

If S^l^ and i^^o known for a model, one-particle, density-matrix, say 7^, then possible 
approximations are the following: 



^co L/mJ(7m=7)- 



(260) 
(261) 



If 7m is equal to, or similar to, the Brueckner-state one-particle density-matrix for an electron 
gas Tg, then, from Eqs. ()25()j) and ()26H1 . we get the following approximation: 



9iiTa=y)- 



(262) 



If periodic boundary conditions are used, the Brueckner orbitals are known to be plane 

u 

waves 2, 



85 



90 



so Tg is known. Because of conservation of momentum, there are 



85 



90 



9l|, but, in general. 



many diagrams that are absent in the S^l''''^[Tg] expansion 
these diagrams appear in the expansion of i^cr'^'W- Since, apparently, all of these excluded- 
diagrams are also excluded in the -E^q"'°'[7] expansion, the above approximation, Eq. ()262j) . 
is probably most appropriate for [i] = III): 



(263) 



For systems with a non-constant external potential, i?^o°'[7] is, in many instances, the 
dominant portion of -Eco'[7) "^^l- ^^at case, we can neglect all terms except the electron gas 
one from Eq. ()258b|) . and using Eq ()263p . we get the following electron-gas approximation: 

ETb]^S^r%](r,=,), (264) 

where, in this approximation, the dependence on v is neglected. 

Eq. ()264|) shares many similarities with the LDA of density-functional theory 
where this approach constructs approximate energy-functionals from expressions derived 
from a uniform electron gas. If periodic boundary conditions are used, the Brueckner orbitals 
(and the Hartree-Fock ones) are known to be plane waves 0, |^, so Tg is known. If, in 
addition to requiring the volume to be infinitely large, the number of particles becomes 
infinite, a uniform electron gas is obtained. In this limiting case, the density of the Brueckner 



reference-state |9), say pug, is identical to the density of the target state both being 
a constant; the correlation energ y o f a uniform electron gas, say ^^cr'^HPug), is a function 
of this density, not a functional jS^ 



85, ^JJ. In the LDA a functional is constructed 



using the function i^cr'^HPug) divided by the number of electrons - the correlation energy 
per particle. An analogous approach may be necessary when constructing the functional 
^riilj9ltu=7); although the one-particle density matrix for an electron gas is not a constant 

ayimi 



|9l|. Furthermore, when eva^ 



85 



uating the diagrams for £^l'"'^{pug), as in the random 



9G 



9l|, the summations over the occupied, plane- wave 



phase approximation (RPA) j54, 
states are replaced by integrals. For an exact treatment of S^l''''^[Tg], this approach cannot be 
used, and, mathematically speaking, this is the difference between i^cr^'f'^g] ^cr'KPug)- 
In using the above electron gas approximation, Eq. (j264|) . or other similar approximations, 
especially for atoms or molecules, one must take into consideration the boundary conditions 
imposed on the wavefunction. For example, if the wavefunction for an electron gas is required 
to vanish at the end points of a cube - and the cube has a finite size - this may lead to 
an electron gas functional £^l''''^[Tg](^rg='r) that explicitly depends on the boundary conditions. 
However, probably in many cases, this explicit dependents can be transformed, or contained 
within Tg, since Tg also depends on the boundary conditions. 

We have, so far, used the external-potential expansions, Eqs. fl258p . to obtain approximate 
functionals. An alternative approach is obtained, if E^^l^['y,v] is known for some system, say 
the helium atom, in the vicinity of some one-particle density-matrix, say the Brueckner one. 
In that case, the following prescription yields an approximate correlation-energy functional: 



where t^^ is the Brueckner, one-particle, density-matrix for the helium atom, and v^^ is the 
external potential for this system, where using Eq. (f75|). we have 

In the limit of p — > r^., necessarily, many terms from Sco[tx,Vx] must vanish, where is 
the external potential associated with the Brueckner one-particle density matrix, Tx- Since, 
apparently, all of these vanishing-diagrams are also excluded in the E^™'°^['j,Vx] expansion, 
most probably, the above approximation, Eq. ()266j) . is most appropriate for (77 = III). 



Furthermore, since the Brueckner density matrix t„c is approximately equal to the Hartree- 
Fock one, say f„e, we can write 



[7,^] ~ ^co[r«.,t;J(,„,=^, (267) 

Using the expansion for £co, given by Eq. ()258a|) . and neglecting all terms except the electron 
gas one, we get 

[7, ^^i^^^^n^H.] (268) 

The above approximation assumes that the terms arising from the helium potential v^^ are 
small; therefore, including them in the following way should yield only a small error: 

^r[7]^^co[r„.,t;J(^„,=,). (269) 

A well known approximation for Ecoi^'nciVyi^ is given by Colle and Salvetti functional 



19|, say i^co[^Ho]; so we have 

[7] [Tec] (270) 

where we have suppressed any mention of fee, since this functional has no explicit dependence 
on the external potential. However, this functional can still possess an implicit dependence 
on Who, since, for example, its four empirical parameters are determined by using data from 
the helium atom. An improved functional, perhaps, can be obtained by appending a poten- 
tial dependence on these parameters and using the approximation given by Eq. ()267|1 . 

We also mention that the correlation potentials v^^ can be treated in a similar manner as 
the correlation-energy functionals E'^J[^, v\, since they also depend on the external potential: 
f^Qft;]; Equations, or approximations, analogous to Eqs. ()258b|) . ()257b|) . ()264|) . and ()265|) . 
are easily defined. However, we now pursue a different approach, permitting correlation 
potentials to be obtained as functional derivatives of variation correlation-energy functionals. 



XI. VARIATIONAL FORMALISM 

We now introduce variational energy-functionals -E'r;[7]- By functional differentiating these 
functionals with respect to the one-particle density-matrix 7, generalized Fock operators are 
obtained. These operators - denoted by C,'^^ - satisfy the same Brillouin-Brueckner and 



commutation relations, Eqs. ()187a|) and ()188|) . as the corresponding non- variational ones, 

T ■ 

Using our trial wavefunctions from Sec. (j^, we can construct variational energy- 
functionals: 

^.M= ^^;^L,^J^^ =i?i[7] + i^^-'[7], (271) 

where the first-order energy Eilj] is given by Eq. p72j) . our variational correlation-energy 
functionals are 



and we have used Eqs. p23p . ()124a|) and p92p . We also define exchange-correlation (xc) 
energy-functionals: 

E^M = mi] - EAl], (273) 

where the exchange energy E^[y] is given by Eq. ()174|1 . 
The exact energy is also given by 

Comparing this Eq. to Eq. (fTUj) gives 

From the variational theorem, the fourth-trial wavefunctions are equal: 

^:v[7] = i?iv[7], (276) 

and we have 

^r[7] = i?r[7]. (277) 
From the variational theorem, we also have 

Er,[l]>S, (278) 

where Eqs. jZH), (071) . (07H1) . (ITfTll . (l7fH|) . and ((221) give the following equalities that appear 
for the Brueckner one-particle density-matrix r: 



s 


= Er,[T] 


= Er,[r], 


(279) 


£co[r] 


= KlV] 


= E^M 


(280) 




= mr] 


= E^^^[r], 


(281) 



where the latter relations use Eqs. ((TTj), and Eqs. (P7H|) and ^17^ indicate that 
the minimization of E^jl'^] occurs at r. 

We now pursue, in a formal way, the minirnization of -^^[7], using an approach that is 
similar to the procedure used by Parr and Yang 2] in their treatment of Hartree-Fock theory, 
where £',,[7] is subject to the constraint that the one-particle density-matrix comes from a 
single-determinantal state 'y(|$)). This condition is imposed by requiring the one-particle 
density- matrix 7 to have a trace equal to the number of electrons AL, and that it is also 
idempotent j49|, |50|| : 



7(X3, X4)(5(X4 - X3) dxs (iX4 = A^-^, 
7(X3,X5)7(X5,X4) dX5 = 7(X3,X4). 



(282) 
(283) 



The normalization constraint, Eq. (j282|) . insures that 7 can be constructed from A'^ orbitals; 
Eq. ()283|1 insures that the density-matrix operator 7 - when acting within the one-particle 
Hilbert space - is a projector into the occupied subspace, as indicated by Eq. ()169|) . where 
7 serves as the kernel of the one-particle density-matrix operator 7, Eq. ()180p . 
Using the above constraints, the variational problem is expressed by 



<^^(7)l. = 0, 



(284) 



where 



7(X3, X4)5(X4 - X3) (iX3 (iX4 - A^-^ 



(ix3 (ix4 q;(x4, X3) 



(285) 



J 7(X3,X5)7(X5,X4) rfXs - 7(X3,X4) 



and where a and /3 are the Lagrangian multipliers. Eq. ()284j) is satisfied when the functional 
derivative of C vanishes: 

5^(7) 



(57(x2,xi) 

where the definition of the functional derivative is 



0, 



<5£(7) 



^^^(7) 
(57(x2,xi 



-57(x2,xi)(ixi dx2. 



(286) 



(287) 



Substituting Eq. ^EB into (jml) . yields 

C|'''(xi,X2) - J rfxg [r(xi,X3)a(x3,X2) + a(xi,X3)r(x3,X2)] + a(xi,X2) - /?5(xi - X2) = 0, 

(288) 



where 



C"(Xi,X2) 



(289) 



57(x2,xi) 

Let the two-body functions, a(xi,X2) and C|'''(xi, X2), serve as kernels of operators, a 
and Cr'''; explicitly, we have 

q;0(xi) = J a(xi,X2)0(x2) (ix2, (290) 
e'0(xi) = y"c'(xi,X2)0(x2)dx2. (291) 
Using this notation, the operator form of Eq. (j288j) is given by 

Cl"' - fd - df + d - = 0, (292) 



where f is defined by Eq. 

The identity operator J, defined by Eq. p79p . or expressed by 

I = f + kr, (293) 

can be used to obtain the following relation: 

a — fa — af = I (a — fa — af) I = krakr — faf, (294) 

which we substitute into Eq. ()292|) : this procedure gives 

C*"' - faf + krOikr - /9 = 0, (295) 

and yields the following requirements: 

krCrT = 0, (296a) 
fCr^r = 0, (296b) 

where we have dropped the 77 superscript since, for (7 = r), all operators are equal within 
these blocks: 

krCr^f = KrCr'^T' = ^rCrT, (297a) 

fC^'/t^ = fCr'^K = fCrK- (297b) 



Eq. ()296ap is yet another representation of the Brillouin-Brueckner condition; comparing 
Eqs. (HHZI) and give 

krCrT = kr^rT, (298a) 

fQrkr- = fPrkr, (298b) 

and it is easily verified that the commutation relation, Eq. (jl88|) . also holds for the variational 
one-body operators Cr^ 

'Cr,r = 0. (299) 
An alternative to the exact Hartree-Fock Eq. p89p is 



C'i^^)= E c;!^'.), (300) 

where the appended rj superscripts appear since the occupied-block matrix-elements 
perhaps, depend on r/. However, we can redefine the variational operators to remove 
this dependence, since Eqs. ()296|) still holds. In any event, we assume that C^''' is independent 
of rj and choose orbitals that diagonalize (^^^, giving a generahzed Hartree-Fock Eq. that is 
an alternative to Eq. ()167|1 : 

C.^[(x)=C#(x). (301) 
Substituting Eq. (l7fT|l into for (r = 7) gives 

C^'(X1,X2) =F^(X1,X2) + Z/-(X1,X2), (302) 

where F^(xi, X2) and ^'^^^(xi, X2) are the kernel of the Fock operator - defined by Eq. (j229p 
- and variational correlation-potentials vlo'- 

F,(xi,X2) = (303) 



^c^'(xi,X2) = , (304) 



57(x2,xi)' 

5 mi] _ 

57(x2,Xi)' 

The operator form of Eq. ()302|) is 

C^' = F^ + z>-. (305) 

Substituting this Eq. into Eq. ()297|) indicates that we have 

KK2r = KK^T = kri'l.f, (306a) 
rK2i^r = ru^okr = fPl^K. (306b) 



Substituting Eqs. ()3U5p and p91|) into ()298p and using the two above definitions, yields 

^ri>;o^ = f^rvl^r, (307a) 
fPl^kr. = fvlX- (307b) 

In order acquire to the kernels of the generalized Fock operators Ct''K^1)^2), given by 
Eq. ()302|) . it is necessary obtain the functional derivatives of Ei[^] and ^^2' [7], as indicated 
by Eq. ()303|) and ()304p . The functional derivative for Ei [7] can be evaluated using Eqs. p72p . 
()173j) . and p74|l . yielding Eq. ()229|1 . The functional derivative of the diagrammatic terms of 
-^co [7] ^Iso be obtained; the details are presented elsewhere js^] • Here we only mention 
that by imposing the same occupied and virtual orbital degeneracies as in Sec. IVlll Bl each 
diagram in the expansion is given by a product of one-particle density-matrices, and can, 
therefore, be differentiated in the same manner as in the treatment of -E'i[7]. After the 
functional derivative is taken, the nondegeneracy of the orbitals can be restored, since the 
entire expansion is invariant to the choice of orbital energies, but this removes the explicit 
dependence on 7 for each term. 

It is easily demonstrated that an analogous external-potential expansion, as given in 
Eq. (|258E|l . also holds: 

m m>n 

Furthermore - by using, Eqs. ()279|) . ()280|) and ()28H) - all approximations presented in Sec.lXl 
are valid when E'^^^ is replaced by £'^2'- For example, consider the electron-gas approximation, 
Eq. (OMl) . and the CoUe-Salvetti functional, Eq (I7fn|l : 

ETb] ^ ^r^N (309) 

^rw ^ ^™[^hJ (310) 

Assuming both approximations are reasonable ones, we can use a linear combination of the 
two: 

^<:,"'[7] ^ ac^S[r„.] + (1 - a,)Slr'['r,\r,=,). (311) 



where Oc is an empirical parameter. This Eq. is an alternative to the B3LYP functional 
which uses analogous correlation-energy functionals: They use the LYP correlation-energy 
functional derived from CoUe-Salvetti one £^l, and a uniform-electron-gas functional, 
derived from the RPA jg^; they set (oc = 8.1). 
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APPENDIX A: VARIATIONAL BRUECKNER-ORBITAL FORMALISMS 



In the exact SCF theory by Lowdin |2^, an orbital variation of an energy functional is 
used to derive the Brillouin-Brueckner condition. By using a slight modification of Kobe's 
formulation tsSj, Lowdin's energy-functional can be written in the following manner: 



I re\ 



(Al) 



where t®, or VqQq, is the reaction operator, which we also define above by the effective 
Hamiltonian Hf^; by definition, this operator satisfies the following relation: 



s = {m!s\<^) = ($|i/(]c,i$),= (1 +x*) 1$), 



(A2) 



and it acts only within the one-dimensional reference-space, with projector P$: {Hfg 
P$/fg^F$); the wave operator and correlation operator x$ are given by Eqs. (0) and 



56 



74 



75 



781, 



79 



931^ 



(I19|) . As in valence-universal multireference perturbation theory 

can be written as a sum of constant, one-, two- and higher-body excitations, where the 
vacuum state is at our disposal. 

By varying the occupied orbitals from the reference state |$), i.e., {i/jq^^}, Lowdin 
and Kobe derived a form of the Brillouin-Brueckner condition, indicating that the above 
functional £',([0,$] - when is held constant - satisfies a variational condition that 
yields the exact energy at the extremum: {S = i?^[0,0]), where the reference state is 
the Brueckner one, |G). A modified version of Lowdin and Kobe's formulation involving a 

n 

one-particle density matrix approach, instead of an orbital one, is easily derived |57| . 



A generalization of Lowdin and Kobe's theory, by Brueckner and Goldman [95", minimizes 
the following functional: i?^ [$,$], given by the right side of Eq. ()A2|) . Such an approach 
has been criticized, since this functional is invariant to the reference state |$), and its a 
constant - the exact energy £\ also, if an approximate Hf^ is used, the energy functional 
does not satisfy the Rayleigh-Ritz principle, as pointed out by Brandow [t^. This method 



leads to the so-called rearrangement potential [95|, |96| that arises from the variation of the 
term involving the reaction operator t*. (For further references regarding the rearrangement 
potential and an historical account of Brueckner orbital theory, see the bibliography notes 
within Kobe's article 

Another variant of the energy functional i?^ [$,$], replaces the correlation operator 
with x% ill the right side of Eq. ()A2|) : this energy functional is simply the non- variational 
one -E'r)['^'], as indicated by Eq. ()130|) . Lowdin and Kobe's exact SCF theory is obtained by 
using Xq and not permitting it to vary, where from Eq. ()132|) we have (xq = Xe)- 



APPENDIX B: PARTITIONING OF SECOND QUANTIZED OPERATORS 



Any second quantized operator, say O, can be partitioned into open (op) and closed (cl) 
portions with respect to a single- or multi-reference space 0, ^3, Ih^- In our case, where 
the reference space is only spanned by a single-determinantal state, the closed portion of O, 
say Och is simply a constant - as in Eqs. (fTTj) and (fT^ - and is given by the fully contracted 
part of O, where the operator is written in normal-ordered form with respect to the reference 



state 1$) j55 



56 



60 



Explicitly, we have 



(Bl) 



where the subscript indicates the zero-body term. 

The open portion of O, say Oop, is usually defined as the remaining portion; it is given 
by the one-, two- and higher-body terms, where, again, the operator is written in normal- 
ordered form. However, for our purposes, we use a more restrictive definition for Oop, and 
define it by the following conditions: 



OopQ# = 0. 



(B2) 
(B3) 



We define the remaining portion, Ore, by the following: 

6 = 0,, + Oop + Ore, (B4) 

so the following identities are satisfied: 

(Oci + 0op)|<f) = 01$), (B5) 

Ore|<f) = 0. (B6) 

The remaining portion O^e has at least one hole or particle annihilation-operator. Oop has 
at least one pair of hole-particle creation-operators and no hole or particle annihilation- 
operators. In terms of diagrams. Ore has at least one external line below the vertex; Oop 
has no lines below the vertex and at least one pair of externals lines above it; Od has no 
external free-lines. 

As an alternative to the above normal-ordered partitioning, we find it convenient to 
partition one-body operators, say h, 

h = J2hij, (B7) 

into the following four components: 

h = he^ + hde + hoc + Kn, (B8a) 

where the excitation (ex), de-excitation (de), occupied (oc), and unoccupied (un) parts are 
given by the following expressions: 

Aex = ^hrwala^,, (B8b) 

wr 

hde = ^/l«;r«I,«r, (B8c) 

wr 

hoc = "^hwxala^, (BSd) 

wx 

Kn = "^hrsaltts. (B8e) 

rs 

and the orbitals are defined with respect to a reference state |$), as indicated by Eqs. ©• 
Note that the open and excited portions are identical: 

^ex = ^op- (B9) 



APPENDIX C: DISCUSSION OF PERTURBATIVE CONVERGENCE WHEN 
USING ORBITAL DEGENERACIES 



Our primary reason for choosing orbital degeneracies within the occupied and virtual 
subspaces is that it provides a means of obtaining an explicit dependence on the one-particle 
density-matrix 7 for the individual perturbative terms (the diagrams) that represent the 
correlation energy £c[l\ and correlation-energy functionals i^^^'W- While this choice greatly 
restricts the zeroth-order Hamiltonian, it can still yield convergent series for ground states, 
as long as the parameter e-y is chosen to be sufficiently large, and negative, so that the 
dominant configurations do not become intruder-states from the presence of small-energy 
denominators or incorrect energy ordering, as in the case for multireference perturbation 



theory 



93 



99 



100 



101 



102 



103[. However, because of incorrect energy-ordering 101 1 



it is probably not possible to generate co nvergent expa nsions for excited states. However, 



71 



93 



104, 



105 



lOq l can be attempted, or the series can 



infinite-order summation methods 
be asymptotically convergent. 

Consider a suitabl e partitioning method, like, for example, M0ller-Plesset 

0, Epstein-Nesbet 



e partitioi 

Hlml, 



112l |. or maximum radius of convergence (i?c) perturbation 



theory 113l . Ill4| . Any of these approaches generate a separate energy-denominator for each 
orthogonal-space state say (Ag = — Eq), where Eq and Eq are the, respective, 
zeroth-order energies. Since small energy-denominators can yield convergence problems, a 
reasonable choice for 2e^ is given by the maximum Ag within the doubly-excited subspace: 

1 



Max Ag, 

2 \q) e m^^)} 



(CI) 



where the search is over all orthogonal-space states {|^)} from within the set of double 
excited-states {I^J^'tr)} that arise from a reasonable set of orbitals, e.g., from the {v^o^ <&,/o} 
and sets. (Some modification is necessary for M0ller~Plesset, since this method 

is restricted to Hartree-Fock orbitals, e.g., choose Hq to be given by Eqs. (fHIIj) . ()38cj) and 
(IHHdl, where (/* = /* = F$).) 

At least for ground states, we anticipate that perturbative expansions using orbital degen- 
eracies can often converge, but, perhaps, at a slow rate. On the other hand, this deficiency is 
partially compensated by a greater computational efficiency, since, for example, costly two- 
electron integral-transformations can be avoided. Furthermore, approximations that involve 



infinite-order summations are often invariant to tlie clioice of Hq - and hence - includ- 
ing tlie coupled cluster method that can be viewed as an infinite-order partial-summation 
method 0,Q- 
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